Introduction to the Renormalization Group 
with Applications to Non-Relativistic 
Quantum Electron Gases 

CIME Lectures, Cetraro 

Vincent Rivasseau 
January 13, 2013 



Abstract 

We review the rigorous work on many Fermions models which lead to 
the first constructions of interacting Fermi liquids in two dimensions, and 
allowed to prove that there are different scaling regimes in two dimensions, 
depending on the shape of the Fermi surface. We also review progress on 
the three dimensional case. 

We start with a pedagogical introduction on quantum field theory and 
perturbative renormalization. Emphasis is then put on using renormal- 
ization around the Fermi surface in a constructive way, in which all orders 
of perturbation theory are summed rigorously. 



1 Introduction to QFT and Renormalization 

1.1 Gaussian Measures 

A finite dimensional centered normalized Gaussian measure d/ic is defined 
through its covariance. Consider a finite dimensional space and a sym- 
metric positive definite N hy N matrix A. The inverse of the matrix A is also 
a definite positive symmetric N hy N matrix C = A^^ called the covariance 
associated to A. The corresponding centered normalized Gaussian measure is 

df,c = (2^)'^/2yd^ e-^*^^^d^X, (1.1) 
so that / djic = 1 • 

To understand c?/ic it is better to know C than A since the moments or 
correlation functions of a Gaussian measure can be expressed simply as sums 



of monomials in C. In fact formula (1.1) perfectly makes sense if C is non 



invertible, and even for C = 0; but the corresponding measure has no density 
with respect to the Lebesgue measure in this case (for C = dfj^c is just Dirac's 



5 function at the origin). The reader famiUar with eg ordinary hnear PDE's 
knows that the essential point is to invert the matrix or the operator, hence 
to know the "Green's function" . But quadratic forms have Hnear equations as 
their variational solutions, so both problems are linked. 

Sine any function can be approximated by polynomials, probability measures 
are characterized by their moments, that is by the integrals they return for each 
polynomial of the integration variables. 

The corresponding theorem which computes the moments of a Gaussian 
measure in terms of the covariance is fundamental in QFT and known there 
under the name of Wick's theorem. It expresses the result as the sum over all 
possible pairings or the variables of a product of covariances between the paired 
variables: 



where G runs over all Wick contractions or pairings of the labels 1, n. Each 
pair i is pictured as a line joining two labels h{t} and e{€} (which we call arbi- 
trarily the "beginning" and "end" of the line) . 

The theory of Gaussian measures and Wick's theorem extends to infinite 
dimensional spaces, in which the covariance C may become a positive kernel 
C{x,y) in a distribution space. We recall that such a kernel is an operator 
acting on functions through C.f = / C{x, y)f{y) dy. The identity operator is 
represented by the Dirac kernel C(x, y) = 5{x — y). But if C is positive definite 
with some regularity it can be also considered the covariance of a truly well 
defined Gaussian measure in some infinite dimensional space of distributions, 
through an extension of Bochner's theorem known as Minlos Theorem [T]. 

1.2 Functional integrals 

In QFT, like in grand-canonical statistical mechanics, particle number is not 
conserved. Cross sections in scattering experiments contain the physical infor- 
mation of the theorjj^ They are the matrix elements of the diffusion matrix S. 
Under suitable conditions they are expressed in terms of the Green functions 
Gm of the theory through so-called "reduction formulae". 

Green functions are time ordered vacuum expectation values of the field 0, 
which is operator valued and acts on the Fock space: 



Here i/ig is the vacuum state and the T-product orders (/)(zi)...0(z7v) according 
to increasing times. 

Consider a Lagrangian field theory, and split the total Lagrangian as the 
sum of a free plus an interacting piece, £ = £o + C,int- The Gell-Mann-Low 

^Correlation functions play this fundamental role in statistical mechanics 
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(1.2) 



Gn{zi, ■■■,zn) =< ipo,T[(l){zi)...(f>{zN)]^o > 



(1.3) 
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formula expresses the Green functions as vacuum expectation values of a similar 
product of free fields with an insertion: 
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V'o > 



(1.4) 



In the functional integral formalism proposed by Feynman [5], the Gell- 
Mann-Low formula is replaced by a functional integral in terms of an (ill-defined) 
"integral over histories" which is formally the product of Lebesgue measures over 
all space time. The corresponding formula is the Feynman-Kac formula: 



Gn{zi, ■■■,Z]y) 
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(1.5) 



The integrand in ( 1.5 1 contains now the full Lagrangian C = Co+Ctnt instead 



of the interacting one. This is interesting to expose symmetries of the theory 
which may not be separate symmetries of the free and interacting Lagrangians, 
for instance gauge symmetries. Perturbation theory and the Feynman rules can 



still be derived as explained in the next subsection. But ( 1.5 ) is also well adapted 



to constrained quantization and to the study of non-perturbative effects. 
For general references on QFT, see OHUl]. 



1.3 Statistical Mechanics and Thermodynamic Quantities 

There is a deep analogy between the Feynman-Kac formula and the formula 
which expresses correlation functions in classical statistical mechanics. 

The partition function of a statistical mechanics grand canonical ensemble 
described by a Hamiltionian H at temperature T and chemical potential fi is 

Za = Tr(e-'^(^-^^)), (1.6) 

where /3 = l/fcT, and the trace may be either a classical integration in phase- 
space or in the quantum case a trace on the relevant Hilbcrt space (Fock space). 

The main problem is to compute the logarithm of the partition function. 
Indeed thermodynamic quantities such as eg the mean energy of the system 



Tr(ffe-'^(-f^-^^)) _ 91ogZA 



< H >T.,= ^-^^ = (1.7) 



the free energy 

F =<H - ^iN >T,^=-\ogZ{,, (1.8) 

the entropy 



1 dF 
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or the heat capacity (at fixed volume) 



fohow from that computation. 

In fact all the detailed information on the equilibrium states is encoded in 
the list of their correlation functions, which are derivatives of the logarithm of 
the partition function with respect to appropriate sources. For instance for a 
lattice Ising model the partition function is 

za= ^"''^"^ (1-11) 

{^x=±l} 

and the correlation functions are 

* ^ {<y.=±i} 

where x labels the discrete sites of the lattice. The sum is over configurations 
{ax = ±1} which associate a "spin" with value +1 or -1 to each such site and 
L{a) contains usually nearest neighbor interactions and possibly a magnetic 
field h: 

x,y nearest neighbors x 



By analytically continuing (1.5) to imaginary time, or Euclidean space, it is 
possible to complete the analogy with (1.12 1, hence to establish a firm contact 
between Euclidean QFT and statistical mechanics O [7l [8] . 



1.4 Schwinger Functions 

This idea also allows to give much better meaning to the path integral, at least 
for a free Bosonic field. Indeed the free Euclidean measure can be defined easily 
as a Gaussian measure, because in Euclidean space Lq is a quadratic form of 
positive typ^ 

The Green functions continued to Euclidean points are called the Schwinger 
functions of the model, and are given by the Euclidean Feynman-Kac formula: 

N 

Sn{z,,...,zn) = Z-' / []0(z,)e-/^(*(-»^-Z?(/., (1.14) 

^However the functional space that supports this measure is not in general a space of 
smooth functions, but rather of distributions. This was already true for functional integrals 
such as those of Brownian motion, which are supported by continuous but not differentiable 
paths. Therefore "functional integrals" in quantum field theory should more appropriately be 
called "distributional integrals". 
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Z= e--''^('*(^»''^ £»<?!). (1.15) 



The simplest interacting field theory is the theory of a one component scalar 
bosonic field </> with quartic interaction ((/)^, which is simpler, is unstable). 
In it is called the (f>\ model. For d = 2, 3 this model is superrenormalizable 
and has been built non perturbatively by constructive field theory (see (SIS]). In 
these dimensions the model is unambiguously related to its perturbation series 
[lOl [TT] through Borel summability fl2|. For d = 4 the model is just renor- 
malizable, and provides the simplest pedagogical introduction to perturbative 
renormalization theory. But because of the Landau ghost or triviality problem 
explained in subsection |1.11[ the model presumably does not exist as a true 
interacting theory at the non perturbative level (see |£, for a discussion of this 
subtle issue). 

Formally the Schwingcr functions of 0^ arc the moments of the measure: 

ie-^-^^'-('"'/2)/0^-(a/2)/(a,0a-0)^^^ (1.16) 



where 



A is the coupling constant, usually assumed positive or complex with pos- 
itive real part; remark the arbitrary but convenient 1/4! factor to take 
into account the symmetry of permutation of all fields at a local vertex. 

TO is the mass, which fixes an energy scale for the theory; 

a is the wave function constant. It can be set to 1 by a rescaling of the 
field. 



• Z is a normalization factor which makes (1.16) a probability measure; 

• D<^ is a formal (mathematically ill-defined) product J| d0(x) of Lebesgue 
measures at every point of . 

The Gaussian part of the measure is 



where Zq is again the normalization factor which makes (1.171 a probability 
measure. 

More precisely if we consider the translation invariant propagator C{x,y) = 
C{x — y) (with slight abuse of notation) , whose Fourier transform is 

g(p)= 2^ 2 ' (1-18) 

we can use Minlos theorem and the general theory of Gaussian processes to de- 
fine dfi{(j)) as the centered Gaussian measure on the Schwartz space of tempered 
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distributions S"(M'') whose covariancc is C. A Gaussian measure is uniquely 
defined by its moments, or the integral of polynomials of fields. Explicitly this 
integral is zero for a monomial of odd degree, and for even n = 2p it is equal to 



/ 



(l){xi)...(t){Xn)dll{(l)) = XI n C'(^bW'^eW), (1-19) 



where the sum runs over all the 2p\\ = {2p — l){2p — 3). ..5. 3.1 Wick pairings W 

of the 2p arguments into the p disjoint pairs £ — (b{i), e{£)). 

Note that since for d>2, C{p) is not integrable, C{x, y) must be understood 
as a distribution. It is therefore convenient to also use regularized kernels, for 
instance 

= 2^ 2 = / e""^' ^"^ (1-20) 

{2'kY P^ + 

whose Fourier transform C^ix,])) is a smooth function and not a distribution: 

C,(x,J/)= re-«™'-(--'')V4a^. (1.21) 

J K ^ 

a~^/'^e~^^~y^'^ is the; Jieat kernel. Therefore this a-representation has also 
an interpretation in terms of Brownian motion: 

C„(a;,j/)= / daexp(— m^a)P(a;, y;a) (1-22) 

J n 

where P{x,y;a) = (47ra)~''/^ exp(— |a; — 2/|^/4a) is the Gaussian probability 
distribution of a Brownian path going from x to j/ in time a. 

Such a regulator k is called an ultraviolet cutoff, and wc have (in the distri- 
bution sense) limK_>.o (^^(a;, y) = C{x,y). Remark that due to the non zero 
mass term, the kernel C^ixjy) decays exponentially at large \x — y\ with rate 
m. For some constant K and > 2 we have: 

\C^ix,y)\ < ii'Kl-'^/2g-m|x-3/|_ (123) 

It is a standard useful construction to build from the Schwinger functions 
the connected Schwinger functions, given by: 



k 

Cn{zu...,zn)= E (-l)'+'n^P*(^ii '•••'%.)' (1-24) 

PlU...UPfc = {l,...,JV}; PinPj=0 i=l 

where the sum is performed over all distinct partitions of {1, A^} into k subsets 
Pi, Pfc, Pi being made of Pi elements called ji, ■■■,jpi- For instance in the (j)'^ 
theory, where all odd Schwinger functions vanish due to the unbroken <p —(f) 
symmetry, the connected 4-point function is simply: 

C4{Z1,...,Z4) = S4{Zi,...,Z4)- S2{Z1,Z2)S2{Z3,Z4) (1.25) 
— S2{zi, Z3)S2{Z2, Z4) — S2{zi,Z4)S2{Z2, Zi). 
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Figure 1: A possible contraction scheme with n = N = 4. 



1.5 Feynman Graphs 

The full interacting measure may now be defined as the multiplication of the 
Gaussian measure d^{<j)) by the interaction factor: 

^e-A-''*''(^)''^d/i(0) (1.26) 
and the Schwinger functions are the normalized moments of this measure: 

S'Ar(2i,...,ZAr) = (j){zi)...(t){zN)di^{(f>). (1-27) 



Expanding the exponential as a power series in the coupling constant A, one 
obtains a formal expansion for the Schwinger functions: 

5^(zi,...,z^) = l£t^ /[/ ^^]>(zi)...0(z^)d/i(0). (1.28) 



Z ^ n\ J 4! 

n=0 

It is now possible to perform explicitly the functional integral of the correspond- 
ing polynomial. The result is at any order n a sum over (4n + A'^ — 1)!! Wick 
contractions schemes W, i.e. over all the ways of pairing together 4ri + A fields 
into 2n + A^/2 pairs. The weight or amplitude of such a scheme W is the spatial 
integral over xi,...,x„ of the integrand Yieew ^(^h{t)^ ■'^■ie(£)) times the factor 
;^(^)". Such amplitudes are functions (in fact distributions) of the external 
positions zi, zn- They may diverge either because they are integrals over all 
of M.^ (no volume cutoff) or because the integrand is typically unbounded due 
to the singularities in the propagator C at coinciding points. 



Labeling the n dummy integration variables in ( 1.28 ) as xi, ...,Xn, we draw a 
line £ for each contraction of two fields. Each position xi , is then associated 
to a four-legged vertex and each external source to a one-legged vertex, as 
shown in Figure [l] 

It is convenient to draw these Wick contractions and to regroup all contrac- 
tions which give rise to the same drawing or graph. There are some subtleties 
about labels. 

Example 1.1 For the normalization at order 1 we have 4 fields, hence 3 Wick 
contractions, which all give the same graph. For the 2 point function at order 
1 we have 6 fields, and 15 Wick contractions which fall into 2 categories with 
weight 3 and 12. 



We have additional observations 



The great advantage of Feynman graphs is that they form a combinatoric 
species in the sense of Joyal ^13] whose logarithm can be computed as the 
species of connected graphs. As we already remarked, the computation of 
this logarithm is the key physical problem. 

However Feynman graphs proliferate, that is their generating functional 
26^0 radius of convergence in A. At the heart of any con- 
structive strategy [H El [HI [HI HH Ull HH], lies the replacement of the 
proliferating species of Feynman graphs by a better one [12] , typically the 
species of forests. The corresponding connected species is the species of 
trees, which does not proliferate. Indeed by Cayley's theorem there are 
only rt"~^ labeled trees on n vertices. This is why constructive expan- 
sions converge while ordinary perturbative expansions dont. Constructive 
theory ultimately may be considered just as repacking Feynman graphs in 
some clever way according to underlying forests [20]. See also the discus- 
sion in subsection 11.161 and below. 

The computation factorizes nicely into the connected components of the 
graphs. These components may or may not have external arguments. In 
the expansion for the normalized functions the vacuum components (i.e. 
those without external arguments) factor out and disappear. Only graphs 
whose connected components all contain external arguments remain. 

If we further search for elementary bricks of the expansion, we can consider 



the connected Schwinger functions like (1.25). In the expansion of these 
functions only the graphs with a single connected component containing 
all external arguments survive. 

1.6 Feynman Rules 

The "Feynman rules" summarize how to compute the amplitude associated to 
a Feynman graph with its correct combinatoric factor. 
We always use the following notations for a graph G: 

• n(G) or simply n is the number of internal vertices of G, or the order of 
the graph. 

• 1{G) or I is the number of internal lines of G, i.e. lines hooked at both 
ends to an internal vertex of G. 



N{G) or N is the number of external vertices of G; it corresponds to the 
order of the Schwinger function one is looking at. When iV = the graph 
is a vacuum graph, otherwise it is called an iV-point graph. 

c{G) or c is the number of connected components of G, 

L(G) or L is the number of independent loops of G. 
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For a regular (f>^ graph, i.e. a graph which has no hne hooked at both ends 
to external vertices, we have the relations: 



1{G) = 2n{G) - N[G) /2, (1.29) 

L{G) = 1{G) ~ n{G) + c(G) = n(G) + 1 - iV(G)/2. (1.30) 

where in the last equality we assume connectedness of G, hence c(G) = 1. 

A subgraph of a graph G is a subset of internal lines of G, together with 
the corresponding attached vertices. Lines in the subset defining F are the 
internal lines of F, and their number is simply 1{F), as before. Similarly all the 
vertices of G hooked to at least one of these internal lines of F are called the 
internal vertices of F and considered to be in F; their number by definition is 
n{F). Finally a good convention is to call external half- line of F every half-line 
of G which is not in F but which is hooked to a vertex of F; it is then the 
number of such external half-lines which we call N{F). With these conventions 



one has for t/)^ subgraphs the same relation (1.291 as for regular graphs. 

To compute the amplitude associated to a (j)* graph, we have to add the 
contributions of the corresponding contraction schemes. This is summarized by 
the "Feynman rules" : 

• To each line £ with end vertices at positions xi and yi, associate a propa- 
gator G{xt,yi). 

• To each internal vertex, associate (— A)/4!. 

• Count all the contraction schemes giving this diagram. The number should 
be of the form {4\)"n\/ S{G) where S{G) is an integer called the symmetry 
factor of the diagram. The 4! represents the permutation of the fields 
hooked to an internal vertex. 

• Multiply all these factors, divide by n! and sum over the position of all 
internal vertices. 

The formula for the bare amplitude of a graph is therefore, as a distribution 
in zi, ....zn: 

/n 
Y[ dxi Y[ C{ 
i=l ieG 

This is the "direct" or "x-space" representation of a Feynman integral. As 
stated above, this integral suffers of possible divergences. But the corresponding 
quantity with both volume cutoff A and ultraviolet cutoff k, namely 

JI^Xj Jl GK(x,,(,),a;j^,,)), (1.32) 

i=l ItzQ 

is well defined. The integrand is indeed bounded and the integration domain A 
is assumed compact. 
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The unnormalized Schwinger functions are therefore formally given by the 
sum over all graphs with the right number of external lines of the corresponding 
Feynman amplitudes: 

ZSn= E ^ gL Ag. (1.33) 

0* graphs G with N(G)=N ^ ' 

Z itself, the normalization, is given by the sum of all vacuum amplitudes: 

(_A)"(G) 



E 



S(G) 

0* graphs G with N {G)=0 ^ ' 



-Ag. (1.34) 



We already remarked that the species of Feynman graphs proliferate at large 
orders. More precisely the total number of 0^ Feynman graphs at order n with 
N external arguments is {An + N)\\. Taking into account Stirling's formula and 
the symmetry factor \/n\ from the exponential we expect perturbation theory 
at large order to behave as i4r"?i! for some constant K . Indeed at order n the 
amplitude of a Feynman graph is a 4n-dimensional integral. It is reasonable 
to expect that in average it should behave as c" for some constant c. But this 



means that one should expect zero radius of convergence for the series (1.33). 
This is not too surprising. Even the one-dimensional integral 

F{g) = / e--'/2-^-'/4'da; (1.35) 



is well-defined only for A > 0. We cannot hope infinite dimensional functional 
integrals of the same kind to behave better than this one dimensional integral. 
In mathematically precise terms, F is not analytic near A = 0, but only Borel 
summable. Borel summability |12j is therefore the best we can hope for the 4>'^ 
theory, and we mentioned that it has indeed been established for the 4)^ theory 
in dimensions 2 and 3 [101 HI] • 

From translation invariance, we do not expect Aq ^ to have a limit as A — ;> oo 
if there are vacuum subgraphs in G. But obviously an amplitude factorizes as 
the product of the amplitudes of its connected components. 

With simple combinatoric verification at the level of contraction schemes we 
can factorize the sum over all vacuum graphs in the expansion of unnormalized 
Schwinger functions, hence get for the normalized functions a formula analog to 



(1.331 



0* graphs G with N{G)=N 
G without any vacuum subgraph 



Now in (1.36) it is possible to pass to the thermodynamic limit (in the sense of 
formal power series) because using the exponential decrease of the propagator, 
each individual graph has a limit at fixed external arguments. There is of course 
no need to divide by the volume for that because each connected component 
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in (1.36) is tied to at least one external source, and they provide the necessary 
breaking of translation invariance. 

Finally one can find the perturbative expansions for the connected Schwinger 
functions and the vertex functions. As expected, the connected Schwinger func- 
tions are given by sums over connected amplitudes: 

0" connected graphs G with N{G)=N ^ ' 

and the vertex functions are the sums of the amputated amplitudes for proper 
graphs, also called one-particle-irreducible. They are the graphs which remain 
connected even after removal of any given internal line. The amputated am- 
plitudes are defined in momentum space by omitting the Fourier transform of 
the propagators of the external lines. It is therefore convenient to write these 
amplitudes in the so-called momentum representation: 

r^(zi,...,ZAr) = siG) ^S(^i.-.^iv), (1.38) 

0-4 proper graphs G with N{G)=N ^ ' 

Ag(zi,...,zjv) ^ (27r)^''^/2 / <^P^-dPNe'^P''''AG{pi,-,PN), (1-39) 
Ag{pu-,Pn)- j n J^.X{^C£<v,l)p,). (1.40) 

£ internal line of G t)GG £ 



Remark in (1.40) the 5 functions which ensure momentum conservation at 
each internal vertex v\ the sum inside is over both internal and external mo- 
menta; each internal line is oriented in an arbitrary way (from 6(£) to e{t)) and 
each external line is oriented towards the inside of the graph. The incidence 
matrix e{v,tj captures in a nice way the information on the internal line^ It 
is 1 if the line i arrives at u, -1 if it starts from v and otherwise. Remark also 
that there is an overall momentum conservation rule 5{pi -I- ... +Pn) hidden in 



(1.40). The drawback of the momentum representation lies in the necessity for 
practical computations to eliminate the i5 functions by a "momentum routing" 
prescription, and there is no canonical choice for that. Although this is rarely 
explicitly explained in the quantum field theory literature, such a choice of a 
momentum routing is equivalent to the choice of a particular spanning tree of 
the graph. 



1.7 Scale Analysis and Renormalization 

In order to analyze the ultraviolet or short distance limit according to the renor- 
malization group method |22j , we can cut the propagator C into slices Ct so that 

•^Strictly speaking this is true only for semi-regular graphs, i.e. graphs without tadpoles, 
i.e. without lines which start and end at the same vertex, see [21j . 



11 



C = X^i^o ^i- This can be done conveniently within the parametric represen- 
tation, since a in this representation roughly corresponds to l/p^- So we can 
define the propagator within a slice as 

(1-41) 

where M is a fixed number, for instance 10, or 2, or e. We can intuitively 
imagine Cj as the piece of the field oscillating with Fourier momenta essentially 
of size M'. In fact it is easy to prove the bound (for d > 2) 

\Ci{x,y)\ < ;i'.M(<^-2)'e-^>-^l (1.42) 

where K is some constant. 

Now the full propagator with ultraviolet cutoff M^, p being a large integer, 
may be viewed as a sum of slices: 

i=0 

Then the basic renormalization group step is made of two main operations: 

• A functional integration 

• The computation of a logarithm 

Indeed decomposing a covariance in a Gaussian process corresponds to a 
decomposition of the field into independent Gaussian random variables each 
distributed with a measure d/Xj of covariance Cj. Let us introduce 

i 

^i = J2'Pr (1-44) 

This is the "low- momentum" field for all frequencies lower than i. The RG 
idea is that starting from scale p and performing p — i steps, one arrives at an 
effective action for the remaining field $i. Then, writing <I>i = + one 
splits the field into a "fiuctuation" field and a "background" field The 
first step, functional integration, is performed solely on the fiuctuation field, so 
it computes 

Z,_i($,_i) = J dfii{<i>i)e-^^^'^^+^^-^K (1.45) 

Then the second step rewrites this quantity as the exponential of an effective 
action, hence simply computes 

5i_i($,_i) = - log[Z,_i($,_i)] (1.46) 

Now Zi-i = e~'^'-i and one can iterate! The flow from the initial bare action 
S = Sp for the full fleld to an effective renormalized action So for the last "slowly 
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Figure 2: A high energy subgraph S seen from lower energies looks quasi- local. 

varying" component 0o of the field is similar to the flow of a dynamical system. 
Its evolution is decomposed into a sequence of discrete steps from Si to Si-i. 

This renormalization group strategy can be best understood on the system 
of Feynman graphs which represent the perturbative expansion of the theory. 
The first step, functional integration over fluctuation fields, means that we have 
to consider subgraphs with all their internal lines in higher slices than any of 
their external lines. The second step, taking the logarithm, means that we have 
to consider only connected such subgraphs. We call such connected subgraphs 
quasi-local. Renormalizability is then a non trivial result that combines locality 
and power counting for these quasi-local subgraphs. 



1.8 Locality, Power Counting 

Locality simply means that quasi-local subgraphs S look local when seen through 
their external lines. Indeed since they are connected and since their internal lines 
have scale say > i, all the internal vertices are roughly at distance Af^*. But 
the external lines have scales < « — 1, which only distinguish details larger than 
Af~(*~^\ Therefore they cannot distinguish the internal vertices of S one from 
the other. Hence quasi-local subgraphs look like "fat dots" when seen through 
their external lines, see Figure [2] Obviously this locality principle is completely 
independent of dimension. 

Power counting is a rough estimate which compares the size of a fat dot 
such as S in Figure [2] with N external legs to the coupling constant that would 
be in front of an exactly local J (j)^ {x)dx interaction term if it were in the 
Lagrangian. To simplify we now assume that the internal scales are all equal 
to i, the external scales are 0(1), and we do not care about constants and so 
on, but only about the dependence in i as i gets large. We must first save 
one internal position such as the barycenter of the fat dot or the position of a 
particular internal vertex to represent the J dx integration in J(j)^{x)dx. Then 
we must integrate over the positions of all internal vertices of the subgraph save 
that one. This brings about a weig ht because since S is connected 

we can use the decay of the internal lines to evaluate these n — 1 i ntegrals. 



Finally we should not forget the prefactor M^^ coming from (1.42), for the 
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I internal lines. Multiplying these two factors and using relation (1.29 )-( 1.30 1 we 
obtain that the "coupling constant" or factor in front of the fat dot is of order 
]^-di{n-i)+2ii2n-N/2) ^ M'^^^\ if wc define the superficial degree of divergence 
of a 0^ connected graph as: 

u;{G) = {d-4)n{G)+d-^^NiG). (1.47) 

So power counting, in contrast with locality, depends on the space-time dimen- 
sion. 

Let us return to the concrete example of Figure [2j A 4-point subgraph made 
of three vertices and four internal lines at a high slice i index. If we suppose 
the four external dashed lines have much lower index, say of order unity, the 
subgraph looks almost local, like a fat dot at this unit scale. We have to save one 
vertex integration for the position of the fat dot. Hence the coupling constant 
of this fat dot is made of two vertex integrations and the four weights of the 
internal lines (in order not to forget these internal line factors we kept internal 
lines apparent as four tadpoles attached to the fat dot in the right of Figure [2]) . 
In dimension 4 this total weight turns out to be independent of the scale. 



1.9 Renormalization, Effective Constants 

At lower scales propagators can branch either through the initial bare coupling 
or through any such fat dot in all possible ways because of the combinatorial 
rules of functional integration. Hence they feel effectively a new coupling which 
is the sum of the bare coupling plus all the fat dot corrections coming from 
higher scales. To compute these new couplings only graphs with w(G) > 0, 
which are called primitively divergent, really matter because their weight does 
not decrease as the gap i increases. 

- If d = 2, we find w(G) = 2 — 2n, so the only primitively divergent graphs 
have n = 1, and = or = 2. The only divergence is due to the "tadpole" 
loop / (p2'^^2 ) which is logarithmically divergent. 

- If d = 3, we find uj{G) = 3 — n — N/2, so the only primitively divergent 
graphs have n < 3, A^ = 0, or n < 2 and N = 2. Such a theory with only a finite 
number of "primitively divergent" subgraphs is called superrenormalizable. 

- If d = 4, uj{G) = 4 — N. Every two point graph is quadratically divergent 
and every four point graph is logarithmically divergent. This is in agreement 
with the superficial degree of these graphs being respectively 2 and 0. The 
couplings that do not decay with i all correspond to terms that were already 
present in the Lagrangian, namely / (j)'^, J (j)^ and /(V(/)).(V(/)^ Hence the 
structure of the Lagrangian resists under change of scale, although the values 
of the coefficients can change. The theory is called just renormalizable. 

^Because the graphs with N = 2 are quadratically divergent we must Taylor expand the 
quasi local fat dots until we get convergent effects. Using parity and rotational symmetry, 
this generates only a logarithmically divergent J(V(/)).(V</)) term beyond the quadratically 
divergent J 0^. Furthermore this term starts only at n = 2 or two loops, because the first 
tadpole graph at A' = 2, n = 1 is exactly local. 



14 



- Finally for c? > 4 we have infinitely many primitively divergent graphs 
with arbitrarily large number of external legs, and the theory is called non- 
renormalizable, because fat dots with N larger than 4 are important and they 
correspond to new couplings generated by the renormalization group which are 
not present in the initial bare Lagrangian. 

To summarize: 

• Locality means that quasi-local subgraphs look local when seen through 
their external lines. It holds in any dimension. 

• Power counting gives the rough size of the new couplings associated to 
these subgraphs as a function of their number N of external legs, of their 
order n and of the dimension of space time d. 

• Renormalizability (in the ultraviolet regime) holds if the structure of the 
Lagrangian resists under change of scale, although the values of the coef- 
ficients or coupling constants may change. For 0^ it occurs if d < 4, with 
d — 4 the most interesting case. 



1.10 The BPHZ Theorem 

The BPHZ theorem is both a brilliant historic piece of mathematical physics 
which gives precise mathematical meaning to the notion of renormalizability, 
using the mathematics of formal power series, but it is also ultimately a dead 
end and a bad way to understand and express renormalization. Let us try to 
explain both statements. 

For the massive Euclidean (/)| theory we could for instance state the following 
normalization conditions on the connected functions in momentum space at zero 
momenta: 

C4(0,0,0,0) = -A,en, (1.48) 
C2(p2 = o) = ^, (1.49) 

HI 

-^CX.=o = -^. (1.50) 
dp^ mj^n 

Usually one puts Uren = 1 by rescaling the field (j). 

Using the inversion theorem on formal power series for any fixed ultraviolet 
cutoff K it is possible to reexpress any formal power series in Xbare with bare 
propagators l/{abarep'^ + iT^lare) f*-"^ ^'^y Schwinger functions as a formal power 
series in A,-e„ with renormalized propagators l/{arenp'^ + ^ren)- The BPHZ 
theorem then states that that formal perturbative formal power series has finite 
coefficients order by order when the ultraviolet cutoff k is lifted. The first proof 
by Hepp (5S] relied on the inductive Bogoliubov's recursion scheme [?]. Then 
a completely explicit expression for the coefficients of the renormalized series 
was written by Zimmermann and many followers |27| . The coefficients of that 
renormalized series can ne written as sums of renormalized Feynman ampli- 
tudes. They are similar to Feynman integrals but with additional subtractions 
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indexed by Zimmermann's forests. Returning to an inductive rather than ex- 
plicit scheme, Polchinski remarked that it is possible to also deduce the BPHZ 
theorem from a renormalization group equation and inductive bounds which 
does not decompose each order of perturbation theory into Feynman graphs 
pS] . This method was clarified and applied by C. Kopper and coworkers, see 
[?]• 

The solution of the difficult "overlapping" divergence problem through Bo- 
goliubov's or Polchinski's recursions and Zimmermann's forests becomes par- 
ticularly clear in the parametric representation using Hepp's sectors. A Hepp 
sector is simply a complete ordering of the a parameters for all the lines of the 
graph. In each sector there is a different classification of forests into packets so 
that each packet gives a finite integral [301 1311 132] • 

But from the physical point of view we cannot conceal the fact that purely 
perturbative renormalization theory is not very satisfying. At least two facts 
hint at a better theory which lies behind: 

• The forest formula seems unnecessarily complicated, with too many terms. 
For instance in any given Hepp sector only one particular packet of forests 
is really necessary to make the renormalized amplitude finite, the one 
which corresponds to the quasi-local divergent subgraphs of that sector. 
The other packets seem useless, a little bit like "junk DNA". They are 
there just because they are necessary for other sectors. This does not look 
optimal. 

• The theory makes renormalized amplitudes finite, but at tremendous cost! 
The size of some of these renormalized amplitudes becomes unreasonably 
large as the size of the graph increases. This phenomenon is called the 
"renormalon problem" . For instance it is easy to check that the renormal- 
ized amplitude (at external momenta) of the graphs P„ with 6 external 
legs and 2 internal vertices in Figure [3] becomes as large as c"n! when 
n — )■ oo. Indeed at large q the renormalized amplitude in Figure [H] 
grows like log \q\. Therefore the chain of n such graphs in Figurejsjbehaves 
as [log IqI]", and the total amplitude of P„ behaves as 



So after renormalization some families of graphs acquire so large values 
that they cannot be resummed! Physically this is just as bad as if infinities 
were still there. 

These two hints are in fact linked. As their name indicates, renormalons are 
due to renormalization. Families of completely convergent graphs such as the 
graphs Qn of Figure [4] are bounded by c", and produce no renormalons. 

Studying more carefully renormalization in the a parametric representation 
one can check that renormalons are solely due to the forests packets that we 
compared to "junk DNA". Renormalons are due to subtractions that are not 




(1.51) 
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Figure 3: A family of graphs P„ producing a renormalon. 




Figure 4: A family of convergent graphs Qn, that do not produce any renor- 
malon. 

necessary to ensure convergence, just like the strange log \q\ growth of Aq^^ at 
large q is solely due to the counterterm in the region where this counterterm is 
not necessary to make the amplitude finite. 

We can therefore conclude that subtractions are not organized in an optimal 
way by the Bogoliubov recursion. What is wrong from a physical point of view 
in the BPHZ theorem is to use the size of the graph as the relevant parameter 
to organize Bogoliubov's induction. It is rather the size of the line momenta 
that should be used to better organize the renormalization subtractions. 

This leads to the point of view advocated in [S]: neither the bare nor the 
renormalized series are optimal. Perturbation should be organized as a power 
series in an infinite set of effective expansions, which are related through the RG 
flow equation. In the end exactly the same contributions are resummed than in 
the bare or in the renormalized series, but they are regrouped in a much better 
way. 

1.11 The Landau ghost and Asymptotic Freedom 

In the case of 0| only the flow of the coupling constants really matters, because 
the flow of TO and of a for different reasons are not very important in the 
ultraviolet limit: 

- the flow of TO is governed at leading order by the tadpole. The bare mass 
mf corresponding to a finite positive physical mass m^g„ is negative and grows 
as AM^* with the slice index i. But since in the i-th slice is also of order 
M^* but without the A, as long as the coupling A remains small it remains much 
larger than mf. Hence the mass term plays no significant role in the higher 
slices. It was remarked in 9 that because there are no overlapping problem 
associated to IPI two point subgraphs, there is in fact no inconvenience to use 
the full renormalized niren all the way from the bare to renormalized scales, 
with subtractions on IPI two point subgraphs independent of their scale. 
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Gj G2 

Figure 5: The 0^ connected graphs with n = 2, N ~ 4. 

- the flow of a is also not very important, fndeed it really starts at two 
loops because the tadpole is exactly local. So this flow is in fact bounded, and 
generates no renormalons. In fact as again remarked in j9j for theories of the 
(/)| type one might as well use the bare value ahare all the way from bare to 
rcnormalized scales and perform no second Taylor subtraction on any IPI two 
point subgraphs. 

But the physics of (pf in the ultraviolet limit really depends of the flow of A. 
By a simple second order computation there are only 2 connected graphs with 
n = 2 and N = 4 pictured in Figure [5] They govern at leading order the flow 
of the coupling constant. 

In the commutative 0| theory the graph Gi does not contribute to the 
coupling constant flow. This can be seen in many ways, for instance after mass 
renormalization the graph Gi vanishes exactly because it contains a tadpole 
which is not quasi-local but exactly local. One can also remark that the graph 
is one particle reducible. In ordinary translation- invariant, hence momentum- 
conserving theories, one-particle-reducible quasi-local graphs never contribute 
significantly to RG flows. Indeed they become very small when the gap i between 
internal and external scales grows. This is because by momentum conservation 
the momentum of any one-particle-reducible line £ has to be the sum of a finite 
set of external momenta on one of its sides. But a finite sum of small momenta 
remains small and this clashes directly with the fact that £ being internal its 
momentum should grow as the gap i grows. Remark that this is no longer be 
true in non commutative vulcanized (^4^, because that theory is not translation 
invariant, and that's why it will ultimately escape the Landau ghost curse. 

So in 04 the flow is intimately linked to the sign of the graph G2 of Figure [sj 
More precisely, we find that at second order the relation between and Ai_i is 

A,_i ~ A, -/3A? (1.52) 

(remember the minus sign in the exponential of the action), where /3 is a 
constant, namely the asymptotic value of J^j j'/ inf (j j')=i I '^^V^j (^i y)^3' i^^ v) 
when z — > 00. Clearly this constant is positive. So for the normal stable 4)\ 



theory, the relation (1.52) inverts into 



A, ~ A,_i+/3A2_i, (1.53) 

so that fixing the rcnormalized coupling seems to lead at finite z to a large, 
diverging bare coupling, incompatible with perturbation theory. This is the 
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Landau ghost problem, which affects both the (f)f theory and electrodynamics. 
Equivalently if one keeps finite as i gets large, Ao = Xren tends to zero and 
the final effective theory is "trivial" which means it is a free theory without 
interaction, in contradiction with the physical observation e.g. of a coupling 
constant of about 1/137 in electrodynamics. 

But in non-Abelian gauge theories an extra minus sign is created by the 
algebra of the Lie brackets. This surprising discovery has deep consequences. 
The flow relation becomes approximately 

A,; ~ A,_i -/3A,A,_i, (1.54) 

with /? > 0, or, dividing by A^A^-i, 

1/A,~ l/A,_i+/3, (1.55) 

with solution ~ ■ ^ more precise computation to third order in fact 

leads to ^ 

l + Ao(/3z + 7bg* + 0(l))" ^^-^^^ 

Such a theory is called asymptotically free (in the ultraviolet limit) because the 
effective coupling tends to with the cutoff for a finite fixed small renormalized 
coupling. Physically the interaction is turned off at small distances. This theory 
is in agreement with scattering experiments which see a collection of almost free 
particles (quarks and gluons) inside the hadrons at very high energy. This was 
the main initial argument to adopt quantum chromodynamics, a non-Abelian 
gauge theory with SU{3) gauge group, as the theory of strong interactions [?]. 

Remark that in such asymptotically free theories which form the backbone 
of today's standard model, the running coupling constants remain bounded be- 
tween far ultraviolet "bare" scales and the lower energy scale where renormalized 
couplings are measured. Ironically the point of view on early renormalization 
theory as a trick to hide the ultraviolet divergences of QFT into infinite un- 
observable bare parameters could not turn out to be more wrong than in the 
standard model. Indeed the bare coupling constants tend to with the ultravi- 
olet cutoff, and what can be farther from infinity than 0? 

Recently it has been shown to all orders of perturbation theory that there 
should be no Landau ghost but an asymptotically safe fixed point for the similar 
RG flow of the non-commutative Grosse-Wulkenhaar model [34l [35l [36l [37] . 
Therefore this model is a kind of "Ising model" for just renormalizable QFT, 
that is a simple model in which one can presumably fully mathematically control 
at last the phenomenon of renormalization in 11 its aspects. 

1.12 Grassmann representations of determinants and Pfaf- 
fians 

Independent Grassmann variables xij---iXn satisfy complete anticommutation 
relations 

XtXj^-XjXz Vi,j (1.57) 
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so that any function of these variables is a polynomial with highest degree one 
in each variable. 

The rules of Grassmann integration are defined by linearity and 

dXi = 0, / XidXi = 1- 



plus the rule that all dx symbols also anticommute between themselves and with 
all x variables. 

The main important facts are 

• Any function of Grassmann variables is a polynomial with highest degree 
one in each variable. 

• PfafRans and determinants can be nicely written as Grassmann integrals; 

The determinant of any n by n matrix can indeed be expressed as a Grass- 
mann Gaussian integral over 2n independent Grassmann variables which it is 
convenient to name as ipi, . . . ,ipn, "^ij ■ ■ • j "^n, although the bars have nothing 
yet at this stage to do with complex conjugation. The formula is 

det d^p,di;,e- ^^^^o^o . (1.53) 

Remember that for ordinary commuting variables and a positive n by n 
Hermitian matrix M 

-y_n**«-^""""'-aSM^ (1,59) 

In short Grassmann Gaussian measures are simpler than ordinary Gaussian 
measures for two main reasons: 

• Grassmann Gaussian measures are associated to any matrix M, there is 
no positivity requirement for M like for ordinary Gaussian measures. 

• Their normalization directly computes the determinant of M, not the in- 
verse (square-root of) the determinant of M . This is essential in many 
areas where factoring out this determinant is desirable; it explains in par- 
ticular the success of Grassmann and supersymmetric functional integrals 
in the study of disordered systems. 



The stubborn reader which remembers the square-root formula in ( 1.1 1 and 



would like to understand the corresponding "reaZ version" of (1.581 is rewarded 
by the beautiful theory of Pfaffians. Clearly commuting Gaussian real inte- 
grals involve symmetric matrices, but Grassmann Gaussian with only n "real" 
integrals must involve n by n antisymmetric matrices. 

The Pfaffian Pf(A) of an antisymmetric matrix A is defined by 

det A = [Vi{A)f . (1.60) 
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and is known to be polynomila in the coefficients of A. This fact is recovered 
easily by writing it as 



Pf(A) = J dxi-dxne 2:,<, x,A,,x, = J d^^^^ dxne 5 i:.,, 



Indeed we have 



detA= / [|dV'^#^e"^•^'^"^' 



(1.61) 



(1.62) 



Performing the change of variables (which a posteriori justifies the complex 
notation) 



(1.63) 



whose Jacobian is i " , the new variables x and w are again independent Grass- 
mann variables. Now a short computation using Aij = —Aji gives 



det A 



(1.64) 



where we used that n — 2p has to be even and that a factor (—1)^ is generated 
when changing Y[ - dxiduJi into ^ - dx i Jli ^'^i- Equation ( 1.64 1 shows why det A 

□ 



1.611. 



is a perfect square and proves 
A useful Lemma is: 

Lemma 1.2 The determinant of a matrix D + A where D is diagonal and A 
antisymmetric has a " quasi- Pfajfian" representation 



det(D + A) ^ / n 



Proof The proof consists in performing the change of variables ( 1.63 1 and can- 
celing carefully the i factors. □ 
There are also normalized Grassmann Gaussian integrals which may be ex- 
pressed formally as 



dfiM 



JYldil^idipie' 

and again are characterized by their two point function or covariance 

ipitpjdpM = Mij. 



(1.66) 



(1.67) 



plus the Grassmann- Wick rule that rt-point functions are expressed as sum over 
Wick contractions with signs. 

For a much more detailed introduction to the rules of Grassmann calculus 
in QFT, we refer to [41] . 
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1.13 Trees, forests and the parametric representation 

Classical evolution can be expanded perturbatively into sums indexed by trees 
whether in quantum field theory the loops of Feynman graphs are essential. 

The hidden trees of the classical system inside QFT can be revealed only 
under scale analysis, since they do not correspond to ordinary spanning trees of 
the graphs, but to the abstract inclusion relations of short range effects (high 
energy quasi local subgraphs) inside larger ones. This point of view has been 
progressively formalized over the years from Bogoliubov to Zimmermann to 
the most recent formalization by D. Kreimer and A. Connes in terms of Hopf 
algebras. 

But ordinary spanning trees of a connected graph also enter in a fascinating 
way in the computation of its amplitude. Since the heat kernel is quadratic it is 
possible to explicitly compute all spatial integrations in a Feynman amplitude. 
One obtains the so-called parametric representation. The result is expressed in 
terms of topological or so-called "Symanzik" polynomials [231 HI] ■ 

The amplitude of an amputated graph G with external momenta p is, up to 
a normalization, in space-time dimension D: 

Aaip) -5{Y.P) I t/c(«)^/^ 11(6- "'daO. (1.68) 
The first and second Symanzik polynomials Ug and Vg are 

where the first sum is over spanning trees T of G and the second sum is over two 
trees T2, i.e. forests separating the graph in exactly two connected components 
E{T2) and F{T2); the corresponding Euclidean invariant (X]ie_E(T2) ^'*)^ 
momentum conservation, also equal to (X]ig_F(T2) ^'»)^- 

The proof of relations ( 1.69a|L69b ) is a special case of the Tree matrix 



Theorem, which we now explain following 
Theorem 1.3 (Tree Matrix Theorem) Let A be an n by n matrix such that 

n 

EA., =0 Vj . (1.70) 

1=1 

Obviously det A — 0. But let A^^ be the matrix A with line 1 and column 1 
deleted. Then 

detA"=^ J]^*-i^' (1-^1) 

T leT 

where the sum runs over all directed trees on {1, n}, directed away from root 
1. 
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This theorem has both a positivity and a democracy aspect: all trees con- 
tribute with positive, equal weights to the determinant. 

Proof of Theorem [rH We use Grassmann variables to write the determinant 
of a matrix with one line and one raw deleted as a Grassmann integral with two 
corresponding sources: 



detA'' = / (di/'dV') (V'iV'i)e 



-ipAip 



The trick is to use (1.70) to write 



(1.72) 



and to obtain 



detA" = J d^dij (V-i^^exp ^ J2 - V-jOV-jj 



dT^dv- (V'lV'i) n (1 - - ^j-)^^ 

by the Pauli exclusion principle. We now expand to get 



(1.73) 



(1.74) 



(1.75) 



yi=(i,j)eg 



where Q is any subset of [n] x [n] , and we used the notation 



(1.76) 



dcf 



(''j)6e 



The tree matrix theorem then follows from the following 

Lemma 1.4 D,g = unless the graph Q is a tree directed away from 1 in which 
case ilg = 1. 

Proof of the lemma : Trivially, if (i, i) belongs to C/, then the integrand of 
^lg contains a factor ip^ — — and therefore Qg vanishes. 

But the crucial observation is that if there is a loop in Q then again fig = 0. 
This is because then the integrand of flj^,iz contains the factor 



^r(fc) - V'r(l) = (^r(fc) " i'rik-l)) + ' ' ' + {i^r(2) " V'r(l)) 



(2) - H^ril 



(1.78) 



Now, upon inserting this telescoping expansion of the factor V't(*:) ~ V-'t(i) ™to 
the integrand of ilj^.n, the latter breaks into a sum of {k — 1) products. For 
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each of these products, there exists an a S TLIklL such that the factor (V'r(Q) ~ 
V'r(Q-i)) appears twice : once with the + sign from the telescopic expansion of 
(V'T(fc) ~ V'r(i))' ^'^d once more with a + (resp. — ) sign if (T(a), r(a — 1)) (resp. 
{T{a — 1), T{a))) belongs to Again, the Pauli exclusion principle entails that 
ng = 0. 

Now every connected component of Q must contain 1, otherwise there is no 
way to saturate the dipi integration. 

This means that Q has to be a directed tree on {1, ...n}. It remains only to 
see that G has to be directed away from 1, which is not too difhcult. □ 
Relations ( 1.69a|[L69b ) follow rather easily from the tree matrix theorem 



and the direct representation of Feynman amplitudes (1.80). 

In [2T] a deeper proof of these relations is given. It relies on the more 
canonical phase-space parametric representation^ which we briefly describe now. 
Let us limit ourselves to "semi-regular" graphs, which have no "tadpoles" that 
it no line starting and ending at the same vertex. These graphs (once their 
lines have been oriented in an arbitrary way) are nicely characterized by their 
incidence matrix, which is a regular 1{G) by n{G) matrix e^^, with 

egy = —1 if line (. exits vertex v 

eiy = +1 if line £ enters vertex v 

eiv = otherwise (l-'''9) 

There are also external momenta p/, / = I,-- - ,N, which we could also also 
orient through a matrix efy. 

The momentum parametric representation then writes 

1{G) n{G)-l 
f,v •' 1=1 11 = 1 

But there is no canonical way to solve for the delta functions, something known 
in physics as the procedure of momentum attribution. So it is better to rewrite 
these amplitudes in the phase-space parametric representation 

.1(G) V-1 



/> 1 r., 1 



e=i v=i 
Integrating over momenta leads to the direct space parametric representation: 

'(G) n(G)-l 



n — d/2~ n 

«-l Ctl v=l 



In [3T] it is shown how the representation (1.80) together with the quasi 



Pfafhan representation of Lemma 1.2 leads to deletion-contraction relations for 
the Symanzik polynomials which allow to compute ( 1.69a|[L69b ) from the theory 
of the universal Tutte polynomial. 
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1.14 BKAR Forest Formula 



Since we want to implement Renormalization Group in a non-perturbative or 
constructive way, we need tools to compute connected functions in a non- 
perturbative way, with the right scaling properties for the convergence radius of 
the expansion. For instance in the just renormalizable case, we need a conver- 
gence radius in the coupling constant which is uniform in the scale index. 

The main such tool is a canonical forest formula \38\ I39j which allows to 
package a perturbative expansion in terms of trees rather than Feynman graphs. 
The advantage was already mentioned several times: the species of trees does 
not proliferate [121 [20] at large orders, in contrast with the species of Feynman 
graphs. 

Consider n points; the set of pairs P„ of such points which has n{n — l)/2 
elements £ — for 1 < i < j < n, and a smooth function / of n{n — l)/2 
variables xg, £ £ Vn- Noting dg for the standard forest formula is 

/(l,...,l)=5][n f dwg]{[Y[dg]f)[xfiM)] (1.80) 

where 

• the sum over T is over forests over the n vertices, including the empty one 

• x'^{\wgiY) is the infimum of the wgi for £' in the unique path from i to j in 
J^, where £ = («, j). If there is no such path, a;|^({wf'}) = by definition. 

• The symmetric n by n matrix X'^dw}) defined by Xf^ — 1 and Xf^ ~ 
a;^({w£/}) for 1 < i < j < n is positive. 

This formula can be viewed as a tool to associate to any pair made of a 
graph G and a spanning forest F (Z G a, unique rational number or weight 
w{G,F) between and 1, called the relative weight of T in G. These weights 
are barycentric or percentage factors, ie for any G 

Y.wiG,F)^l (1.81) 

FCG 

The numbers w{G, F) are multiplicative over disjoint unions Hence it is 
enough to give the formula for (G, F) only when G is connected and _F = T is a 
spanning tree in i f\ The definition of these weights is 

Definition 1 ^ 

wiG,T) = l[ f l[dwgllxJi{w}) (1.82) 

where x^i^wY) is again the infimum over the Wgi parameters over the lines £' 
forming the unique path in T joining the ends of £. 

^And also over vertex joints of graphs, just as in the universality theorem for the Tutte 
polynomial 

•^It is enough in faet to compute such weights for 1-particle irreducible and 1-vertex- 
irreducible graphs, then multiply them in the appropriate way for the general case. 
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Consider the expansion in terms of Feynman amplitudes of a connected quan- 
tity S. The most naive way to reorder Feynman perturbation theory according 
to trees rather than graphs is to insert for each graph the relation ( 1.81 ) 



5 = ^Ag = 5]5]«;(G,T)^g 



(1.83) 



G G TCG 



and exchange the order of the sums over S and T. Hence it writes 



(1.84) 



T GDT 



This rearranges the Feynman expansion according to trees, but each tree 
has the same number of vertices as the initial graph. Hence it reshuffles the 
various terms of a given, fixed order of perturbation theory. Remark that if the 
initial graphs have say degree 4 at each vertex, only trees with degree less than 
or equal to 4 occur in the rearranged tree expansion. 

For Fermionic theories this is typically sufficient and one has for small enough 
coupling 



because Fermionic graphs essentially mostly compensate each other at a fixed 
order by Pauli's principle; mathematically this is because these graphs form a 
determinant and the size of a determinant is much less than what its permuta- 
tion expansion usually suggests. This is well known [46l [62l |47] . 

But this repacking fails completely for Bosonic theories, because the only 
compensations there occur between graphs of different orders. Hence if we 
perform this naive reshuffling, eg on the 0g theory we would still have 



Recently a new expansion called the Loop Vertex Expansion has been found 
[16] which overcomes this difficulty by exchanging the role of vertices and prop- 
agators before applying the forest formula. It can also be seen as a combination 
of the forest formula with the so-called intermediate field representation, which 
expands into essentially square roots of a stable Bosonic interaction. We refer 
the reader to [121 [HI [THl dQj i20j but wont review this expansion here, since from 
now on we are mostly going to deal with Fcrmions. 

1.15 Gram and Hadamard Bounds 

These two bounds on a determinant are often confused! 

The Gram bound applies to a matrix A = aij whose entries are scalar 
product. This means we suppose that there exists some Hilbert space H and 
2n vectors fi, i = 1, ■ ■ ■ ,n, gj, j ~ 1, ■ ■ ■ ,n with 




(1.85) 



T 




(1.86) 



T 



aij —< fn 9j >H ■ 



(1.87) 
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The Gram bound states 



idetAi<nii/j//nii5.-ii^- 



(1- 



Of course any matrix A can always be written of the Gram type, eg with 
H = M", fi = (ai^k) and gj = Sj^k, fc = 1, • • • ,n, or conversely. Hence there 
are two corresponding asymmetric Hadamard bounds, one for rows and one for 
columns: 



idetAi <n. ^ 

i=l \ j= 



E4' 



Idet^l < Yl 



(1- 



(1.90) 



and also a symmetric Hadamard bound involving the supremum of the matrix 
elements: 

|detA| < n"/2(sup|a,j|)". (1.91) 

hi 

Remark that the symmetric Hadamard bound means that a determinant of 
a large matrix is always much smaller than what its permutation expansion plus 
naive bounds would suggest, which is the "stupid bound" 



jdet^l < n!(sup |a,y|)". 



(1.92) 



This difference in constructive theory is essential. Indeed for Fermionic theories 
with bounded propagators and a quartic interaction, the matrix A at n-th order 
of perturbation is a 2n x 2n matrix, with propagators as matrix elements, and 



there is a 1/n! symmetry factor. Hence the bound (1.92) would lead to believe 



that the radius of convergence of the partition function is 0, like in the Bosonic 
case. But the Hadamard bound (1.91 1 proves that it is at least finite. Moreover 
usually it is possible to write the propagators as scalar products in 



of 



functions which also have bounded norm^ In that case the Gram bound 
( 1.88 ) shows that the partition function is in fact an entire function, as it shows 
no factorial dependence at all as n -> oo! 



1.16 Single Scale Constructive Theory for a Toy Model 

Consider a just-renormalizable QFT theory. The key problem is to compute 
connected quantities with an expansion which converges for a small coupling 
constant, with a propagator limited to a single renormalization group scale, 
uniformly in the slice index. In the simple Fermionic case, this can be done 
through applying first the BKAR formula then checking convergence through a 
Gram bound. 

'^This is usually easily done by taking some kind of "square roots" in momentum space. 
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To discuss the type of Fermionic models met in condensed matter it is ap- 
propriate to consider first a toy Fermionic d-dimensional QFT model. It is made 
of a single an infrared slice with N colors. Suppose the propagator is diagonal 
in color space and satisfies the bound 

f.r-dj/2 

\Gj,ab{x,v)\ < Sab—^e-''''^^-y^. (1.93) 



We say that the interaction is of the vector type (or Gross-Neveu type) if it 
is of the form 



JV JV 



a=l 6=1 

where A is the coupling constant. 
We claim that 

Lemma 1.5 The perturbation theory for the connected functions of this single 
slice model has a radius of convergence in X which is uniform in j and N . 

To prove this lemma, expands the partition function Z{A) through the forest 
formula, and take the logarithm to obtain a tree formula for the pressure 

p= lim ^logZ(A). (1.95) 

A-i-oo |A| 

This is completely straightforward, the only difficulty being notational. Using 
the notations of [JT] 

p = lim dtici^,^P)e^-^^^'^^^) (1.96) 

A^oo |A| V 



dwi 



n=0 ai,...,a„,6i,...,6„ = l T O leT^ 

/ dxi...dxnS{xi = 0)T\ (C{xin),Xin))Si) X det {Cai3)aeA,0eB 

The sum over the a.^'s and hiS are over the colors of the fields and antifields 
of the vertices obtained by expanding the interaction and of the form: 

i^ai {Xi)lpai {Xi)^bi {Xi)'ll^b, (Xi) (1-97) 

with I < i < n. The sum over T is over all trees which connect together the n 
vertices at a:i, . . . , x„. The sum over D, is over the compatible ways of realizing 
the bonds I — S T as contractions of a and ip between the vertices i and 

j (compatible means that we do not contract twice the same field or antifield). 
e(T, ri) is a sign which is not important for the bound (see |47] for its explicit 
computation). For any I T , i{l) ^ {1, ■ ■ ■ , n} labels the vertex where the field. 
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contracted by the procedure f2 concerning the hnk I, was chosen. Likewise 
is the label for the vertex containing the contracted antifield. Si is 1 if the colors 
(among ai, . . . , a„, &i, . . . , 5„) of the field and antifield contracted by I are the 
same and else is 0. Finally the matrix (Cq^)q_^ of the remaining "loop lines" is 
defined in the following manner. 

The row indices a label the 2n fields produced by the n vertices, so that 
a = (i, cr) with I < i < n and a takes two values 1 or 2 to indicate whether the 



field is the second or the fourth factor in (1.97) respectively. 

The column indices /3 label in the same way the 2n antifields, so that f3 = 
(j, t) with I < j < n and t = 1 or 2 according to whether the antifield is 



the first or the third factor in (1.971 respectively. The a's and f3's are ordered 
lexicographically. We denote by c{i,a) the color of the field labeled by (j,cr) 
that is fli, a a — 1, and 5^ if cr = 2. We introduce the similar notation c{j,T) 
for the color of an antifield. Now 

C{i,a){j,T) wjj^^ {w)C{x,,Xj)S{c{i,a),c{j,T)). (1.98) 

Finally each time a field {i,a) is contracted by the corresponding row is 
deleted from the 2n x 2n matrix (Ca/3). Likewise, for any contracted antifield 
the corresponding column is erased. A and B denote respectively the set of 
remaining rows and the set of remaining columns. The minor determinant 
featuring as detremaining in (1.97) is now det{bap)aeA,peB which is (n + 1) x 



(n + 1). Indeed for each of the n — 1 links of T, a row and a column have been 
erased. 

Suppose we have written 

Cjixk,y„i) < .fj,k,gj,m >L2 (1.99) 

(this is realized through fj^k = fjixk,-) and gj^m = 9j{-,Vm) if fj{p).gj{p) = 
Cj{p)). and that the power counting is conserved by taking square- roots, namely 

that the norms of gj and fj scale as ^^i^^"^ ■ 
Then applying the Gram inequality 

I det[Cj^Q/j()jj.Qj^aining 

anti— fields k fields m 

will lead to the proof that the radius of convergence of the pressure is uniform 
in J and TV for this toy model. 

Indeed the following lemma shows that the presence of the weakening factors 
w does not change the outcome of the Gram bound. 

Lemma Let A = {aap)a,p be a Gram matrix: a^p —< /q,5^ > Jor some inner 

product <.,.>. Suppose each of the indices a and (3 is 0/ the form (i, a) where 
the first index i, 1 < i < n, has the same range as the indices of the positive 



matrix (wf^'^^ {'w))ij , and a runs through some other index set E. 



Let C = {Cap)a,p be the matrix with entries C(i^a)(j.T) — wfj'^^ (w). < 
f(i,a)i g{j,T) > o,n-d let {Cap)aeA.0eB be some square matrix extracted from C, 
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then for any w we have the Gram inequality: 

|det(C„^)„eA,/3eB| < J] H/^H 11 WdpW (1-101) 

Proof Indeed we can take the symmetric square root v of the positive matrix 
^T,BK that w'ij^^ ~ X]fc=i ''^jfe^fej- us denote the components of the 
vectors / and in an orthonormal basis for the scalar product < ., . > with 
q elements, by /"^.-j and g™j ^.y 1 < m < q. (Indeed even if the initial Hilbert 
space is infinite dimensional, the problem is obviously restricted to the finite 
dimensional subspace generated by the finite set of vectors / and g). We then 
define the tensorized vectors -F'(i,CT) and G^j^r) with components F^^™-) = VikfJ^ „^ 
and G^™.-) = "^ikQi] t) where 1 < k < n and 1 < m < q. Now considering the 
tensor scalar product < ., . >t we have 

n q 
k = l m=l 

By Gram's inequality using the < ., . >t scalar product we get 

\det{C^f,)^^A,peB)\ < n II/"IIt R H^/sHt (1-103) 

but 

n q n q 

A; — 1 m— 1 k—1 m— 1 

= ^.»E (/".))' = ll/(.,.)ll', (1-104) 

m— 1 

since Wa = 1 for any «, 1 < i < n. □ 
Let us apply the Gram bound and this last Lemma to bound the determinant 



in (1.97) 



There is a factor M ''■'/^ per line, or M '^^/^ per field ie entry of the loop 
determinant. This gives a factor M~'^^ per vertex. 

There is a factor Af per vertex spatial integration (minus one) 



Hence the A radius of convergence is uniform in j. 

• There is a factor N^^^"^ per line, or N^^^'^ per field ie entry of the loop 
determinant. This gives a factor N^^ per vertex 

• There is a factor N per vertex (plus an extra one) 
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This leads to a bound in N.M"^^ [cA]" for the n-th order of the pressure, hence 
to a radius of convergence at least 1/c. As expected, this is a bound uniform in 
the slice index j. Hence the A radius of convergence is uniform in N . 

The last item, namely the factor A^"+i for the colors sums is the only one 
not obvious to prove. Indeed don't know all the graph, but only a spanning 
tree. We need to organize the sum over the colors from the leaves to the root of 
the tree. In this way the pay a factor N at each leaf to know the color index 
which does not go towards the root, then prune the leaf and iterate. The last 
vertex (the root) is the only special one as it costs two N factors. 

Let us remark that to treat the corresponding toy model in the Bosonic case 
the standard constructive method would be to perform a cluster expansion with 
respect to a lattice of cubes, then a Mayer expansion which further removed the 
remaining hardcore constraints with respect to the cubes [S]. Both expansions 
needed to use the forest formula. This was simplified by the invention of the 
Loop vertex expansion, in which cubes, cluster and Mayer expansions are no 
longer needed. In addition the Loop vertex model leads to uniform bounds also 
for matrix toy models, a result which cannot be obtained up to now with other 
methods [TB] . 

2 Interacting Fermions in Two Dimensions 
2.1 Introduction 

One of the main achievements in renormalization theory has been the extension 
of the renormalization group of Wilson (which analyzes long-range behavior gov- 
erned by simple scaling around the point singularity p = in momentum space) 
to long-range behavior governed by extended singularities [5Tl[52l[55. This very 
natural and general idea is susceptible of many applications in various domains, 
including condensed matter (reviewed here, in which the extended singularity is 
the Fermi surface) but also other ones such as diffusion in Minkowski space (in 
which the extended singularity is the mass shell) . In this section we will discuss 
interacting Fermions models such as those describing the conduction electrons 
in a metal. 

The key features which differentiate electrons in condensed matter from Eu- 
clidean field theory, and make the subject in a way mathematically richer, is 
that space-time rotation invariance is broken, and that particle density is finite. 
This finite density of particles creates the Fermi sea: particles fill states up to 
an energy level called the Fermi surface. 

The field theory formalism is the best tool to isolate fundamental issues such 
as the existence of non-perturbative effects. In this formalism the usual Hamil- 
tonian point of view with operators creating electrons or holes is superseded by 
the more synthetic point of view of anticommuting Fermionic fields with two 
spin indices and arguments in d -f 1 dimensional space-time,. Beware however 
of the QFT-convention to always call dimension the dimension of space-time, 
whether from now on we have to stick to the usual condensed matter conven- 
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tion which is to always call dimension the dimension of space only. So one 
dimensional interacting Fermions correspond at zero temperature to a two di- 
mensional QFT, two dimensional Fermions correspond at zero temperature to 
a three dimensional QFT and so on. 

After the discovery of high temperature superconductivity, a key question 
emerged. Do interacting Fermions in 2 dimensions (above their low-temperature 
phase) resemble more three dimensional Fermions, i.e. the Fermi liquid, or one 
dimensional Fermions, i.e. the Luttinger liquid? The short answer to this con- 
troversial question is that it was solved rigorously by mathematical physics and 
that the answer depends on the shape of the Fermi surface. Interacting Fermions 
with a round Fermi surface behave more like three dimensional Fermi liquids, 
whether interacting Fermions with the square Fermi surface of the Hubbard 
model at half-filling behave more like a one-dimensional Luttinger liquid. 

This statement has been now proved in full mathematical rigor, beyond 
perturbation theory, in the series of works [SU [551 [Ml [53 [53 [53 HO] • 

The existence of usual 2D interacting Fermi liquids was established in [5^155] 
using the mathematically precise criterion of Salmhofer (48, in the case of a 
temperature infra-red regulator. Using a magnetic field regulator that breaks 
parity invariance of the Fermi surface it was also established in the initial sense 
of a discontinuity at the Fermi surface in the series of papers |49j . 



2.2 The Models: Ja, ^3, H2... 

We consider a gas of Fermions in thermal equilibrium at temperature T, with 
coupling constant A. The free propagator for this model is 

Cabik) = Sab ., ^ ... (2.1) 

iKo — e(k) 

with k being the d-dimensional momentum, e(k) — e(k) — /i, e(k) being the 
kinetic energy and fi the chemical potential. The a, 5 G {^, 4,} index is for spin 
hence can take two values (remember spin is treated non-relativistically). 

At finite temperature, since Fermionic fields have to satisfy antiperiodic 



boundary conditions, the component fco in (2.1) can take only discrete values 
(called the Matsubara frequencies) : so the integral over fcg is really a discrete 
sum over n. 

These Matsubara frequencies are: 

fco = ^^^^ , n e Z (2.2) 

where /3 = (kT)^^. For any n we have fco 7^ 0, so that the denominator in C(fc) 
can never be 0. This is why the temperature provides a natural infrared cut-off. 

We can think of fco as the Fourier dual to an imaginary Euclidean-time con- 
tinuous variable taking values in a circle, with length proportional to inverse 
temperature /3. When T — > 0+, (which means f3 — ?> -l-oo), fco becomes a con- 
tinuous variable, the corresponding discrete sum becomes an integral, and the 
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corresponding propagator Cq{x) becomes singular on the Fermi surface defined 
by fco = and e(k) = 0. 

This Fermi surface depends on the kinetic energy e(k) of the modeL 
For rotation invariant models, e(k) = k^/2m where m is some effective 
or "dressed" electron mass. In this case the energy is invariant under spatial 
rotations and the Fermi surface is simply a circle in two dimensions and a sphere 
in three dimensions, with radius ^/2mfi. This jellium isotropic propagator is 
realistic in the limit of weak electron densities. We call this propagator the 
jellium propagator. We always consider this model, the most natural one in the 
continuum, together with an ultraviolet cutoff (which it is natural to also take 
rotation invariant 

Another model considered extensively is the half-filled 2d Hubbard model, 
nicknamed i?2- In this model the position variable x lives on the lattice Z^, and 
e(k) — cos fci + cos ^2 so that at /i = the Fermi surface is a square of side size 
-\/27r, joining the points (tt, 0), (0, tt) in the first Brillouin zone. 

This propagator is called the Hubbard propagator. 



2.3 Interaction, Locality 

The physical interaction between conduction electrons in a solid could be very 
complicated; the naive Coulomb interaction is in fact subject to heavy screening, 
the main effective interaction being due to lattice phonons exchange and other 
effects. But we are interested in long-range physics, Hence we should use a 
quasi-local action which decays It is a bit counterintuitive but in fact perfectly 
reasonable for a mathematical idealization to use in fact a fully local interaction. 
This should capture all essential mathematical difficulties of the corresponding 
renormalization group. 

But there is a unique exactly local such interaction, namely 

Sv = X f d^+' ( V Mx)Mx)) , (2.3) 

where V := [— /?, /3[xT/' and V' is an auxiliary volume cutoff in two dimensional 
space, that will be sent to infinity in the thermodynamic limit. Indeed any local 
polynomial of higher degree is zero since Fermionic fields anticommute. Remark 
it is of the same form than (??), with spin playing the role of color. 

Hence from the mathematical point of view, in contrast with the propagator, 
the interesting condensed matter interaction is essentially unique. 



The models with jellium propagator and such an interaction ( 2.3 1 are respec 



tively nicknamed J2 and J3 in dimensions 2 and 3. The model with Hubbard 



propagator and interaction (2.3 1 is nicknamed i?2- 

It is possible to interpolate continuously between H2 and J2 by varying the 
filling factor of the Hubbard model. Lattice models with next-nearest neighbor 



° The question of whether and how to remove that ultraviolet cutoff has been discussed 
extensively in the literature, but we consider it as unphysical for a non-relativistic model of 
condensed matter, which is certainly an effective theory at best. 
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hopping are also interesting, as they are reahy the ones used to model the high 
Tc superconducting phase in cuprates, but we shall not consider them here for 
simplicity. 

The basic new feature which changes dramatically the power counting of 
the theory is that the singularity of the jellium propagator is of codimension 
2 in the d + 1 dimensional space-time. Instead of changing with dimension, 
like in ordinary field theory, perturbative power counting is now independent 
of the dimension, and is the one of a just renormalizable theory. Indeed in 
a graph with 4 external legs, there are n vertices, 2n — 2 internal lines and 
L — n — 1 independent loops. Each independent loop momentum gives rise to 
two transverse variables, for instance /cq and |k| in the jellium case, and to d — 1 
inessential bounded angular variables. Hence the 2L = 2{n — 1) dimensions of 
integration for the loop momenta exactly balance the 2n — 2 singularities of the 
internal propagators, as is the case in a just renormalizable theory. 

In one spatial dimension, hence two space-time dimensions, the Fermi surface 
reduces to two points, and there is also no proper BCS theory since there is 
no continuous symmetry breaking in two dimensions (by the "Mermin- Wagner 
theorem"). Nevertheless the many Fermion system in one spatial dimension 
gives rise to an interesting non-trivial behavior, called the Luttinger liquid |43j . 

The marvel is that although the renormalization group now pinches a non 
trivial extended singularity, the locality principle still works. The leading part 
of the four point quasi-local graphs is still local, like the interaction of the 
initial theory. This is very surprising since quasi-local graphs have internal 
lines which simply carry excitations farther from the Fermi surface than their 
external ones. Momenta close to the Fermi surface, when moved to a barycenter 
of the graph, should react through a non trivial phase factor. But the miracle 
is that the only divergent part of the main one-loop contribution, the "bubble 
graph" comes when the combination of the two external legs at each end carries 
approximately zero total momentum. This is because only then can the two 
inner bubble propagators both range over the full Fermi sphere. This special 
configuration, with proper spin and arrows to ensure no oscillations occur, is 
really the Cooper pairl 

Then these two external legs with approximately zero total momentum can 
be moved together like a single low momentum leg in an ordinary Wilsonian 
renormalization. This is in essence why renormalization still works in condensed 
matter and beyond the two point function renormalization (Fermi radius renor- 
malization) only changes the value of the coupling constant for the local inter- 
action (2.3). Of course if contributions beyond one loop are taken into account, 
the story becomes more complicated and the renormalization group flow can in 
fact involve infinitely many coupling constants '64' . 

A still more surprising case where locality works in a new form, called Moy- 
ality, is the case of non commutative field theory on Moyal space [31]. There 
again the divergent subgraphs are exactly the only ones that can be renormal- 
ized through counterterms of the initial form of the Lagrangian [5U] . This lead 
us to hope that another still generalized form of locality might hold in quantum 
gravity. 
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2.4 A Brief Review of Rigorous Results 



What did the programs of rigorous mathematical study of interacting Fermi 
systems accomphsh until now? Recall that in dimension 1 there is neither su- 
perconductivity nor extended Fermi surface, and Fermion systems have been 
proved to exhibit Luttinger liquid behavior [iS]. The initial goal of the studies 
in two or three dimensions was to understand the low temperature phase of 
these systems, and in particular to build a rigorous constructive BCS theory of 
superconductivity. The mechanism for the formation of Cooper pairs and the 
main technical tool to use (namely the corresponding 1/N expansion, where N 
is the number of sectors which proliferate near the Fermi surface at low temper- 
atures) have been identified [62| . But the goal of building a completely rigorous 
BCS theory ab initio remains elusive because of the technicalities involved with 
the constructive control of continuous symmetry breaking. So the initial goal 
was replaced by a more modest one, still important in view of the controver- 
sies over the nature of two dimensional Fermi liquids [61 , namely the rigorous 
control of what occurs before pair formation. 

As is well known, sufficiently high magnetic field or temperature are the 
two different ways to break the Cooper pairs and prevent superconductivity. 
Accordingly two approaches were devised for the construction of "Fermi liquids. 
One is based on the use of non-parity invariant Fermi surfaces to prevent pair 
formation. These surfaces occur physically when generic magnetic fields are 
applied to two dimensional Fermi systems. In the large series of papers [45] . 
the construction of two dimensional Fermi liquids for a wide class of non-parity 
invariant Fermi surfaces has been completed in great detail by Feldman, Knorrer 
and Trubowitz. These papers establish Fermi liquid behavior in the traditional 
sense of physics textbooks, namely as a jump of the density of states at the 
Fermi surface at zero temperature, but they do not apply to the simplest Fermi 
surfaces, such as circles or squares, which are parity invariant. 

The other approach is based on Salmhofer's criterion, in which temperature 
is the cutoff which prevents pair formation. The corresponding program studies 
whether given models satisfy Salmhofer's criterion or not. The study of each 
model has been divided into two main steps of roughly equal difficulty, the 
control of convergent contributions and the renormalization of the two point 
functions. In dimension two the corresponding analysis has been completed for 
J2, a Fermi liquid in the sense of Salmhofer, and for H2 which is not, and is a 
Luttinger liquid with logarithmic corrections, according to [54l l55t [56 l [57 | [58] . 

Similar results similar have been also obtained for more general convex curves 
not necessarily rotation invariant such as those of the Hubbard model at low 
filling, where the Fermi surface becomes more and more circular, including an 
improved treatment of the four point functions leading to better constants [E9\ 
[5U] . Therefore as the filling factor of the Hubbard model is moved from half- 
filling to low filling, we conclude that there must be a crossover from Luttinger 
liquid behavior to Fermi liquid behavior. This sheds light on the controversy |61) 
over the Luttinger or Fermi nature of two-dimensional many-Fermion systems 
above their critical temperature. 
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2.5 Multiscale Analysis, Anguleir Sectors 



For any two-dimensional model built until now in the constructive sense, the 
strategy is the same. It is based on some kind of multiscale expansion, which 
keeps a largo fraction of the theory in unexpanded determinants. The global 
bound on these determinant (using determinant inequalities such as Gram in- 
equality) is much better than if the determinant was expanded into Feynman 
graphs which would then be bounded one by one, and the bounds summed. The 
bound obtained in this way would simply diverge at large order (i.e. not prove 
any analyticity at all in the coupling constant) simply because there are too 
many Feynman graphs at large order. But the divergence of a bound does not 
mean the divergence of the true quantity if the bound is bad. Constructive anal- 
ysis, which keeps loops unexpanded is the correct way to obtain better bounds, 
which do prove that the true series in fact does not diverge, i.e. has a finite 
convergence radius in the coupling constant. This radius however shrink when 
the temperature goes to 0, and a good constructive analysis should establish 
the correct shrinking rate, which is logarithmic. This is where multiscale rather 
than single scale constructive analysis becomes necessary. 

The basic idea of the multiscale analysis is to slice the propagator according 
to the size of its denominator so that the slice with index j corresponds to 
\iko + e(k)| ~ M~\ where M is some fixed constant. 

This multiscale analysis is supplemented within each scale by an angular 
"sector analysis. The number of sectors should be kept as small as possible, so 
each sector should be as large as possible in the directions tangent to the Fermi 
surface in three dimensions, or to the Fermi curve in two dimensions. What 
limits however the size of these sectors is the curvature of the surface, so that 
stationary phase method could still relate the spatial decay of a propagator 
within a sector to its dual size in momentum space. In the case of a circle, the 
number of sectors at distance of the singularity grows therefore at least 

like M-'/^, hence like a power of T. However for the half-filled Hubbard model, 
since the curvature is "concentrated at the corners the number of sectors grows 
only like | logTj. In one dimension there are really only two sectors since the 
Fermi singularity is made of two points. A logarithm is closer to a constant 
than to a power; this observation is the main reason for which the half-filled 
Hubbard model is closer to the one-dimensional Luttinger liquid than to the 
three dimensional Fermi liquid. 

Momentum conservation rulcis for sectors which meet at a given vertex in 
general are needed to fix the correct power counting of the subgraphs of the 
model. In the Hubbard case at half filling, these rules are needed only to fix 
the correct logarithmic power counting, since the growth of sectors near the 
singularity is only logarithmic. In both cases the net effect in two dimensions of 
these conservation rules is to roughly identify two pairs of conserved sectors at 
any vertex, so that in each slice the model resembles an A''- component vector 
model, where N is the number of sectors in the slice. 

The multiscale renormalization group analysis of the model then consists 
essentially in selecting, for any graph, a tree which is a subtree in each of the 
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quasi-local connected components of the graph accord- ing to the momentum 
slicing. These connected components are those for which all internal lines are 
farther from the Fermi surface than all external lines. The selection of this tree 
can be performed in a constructive manner, keeping the remaining loop fields 
in a determinant. The combinatoric difficulty related to the fact that a graph 
contains many trees is tackled by the forest formula. 

Once the scale analysis has been performed, a partial expansion of the loop 
determinant can detect all the dangerous two and four point functions which 
require renormalization. A key point is that this expansion can be done without 
destroying the Gram bound, and the corresponding sum is not too big (this 
means its cardinal remains bounded by (where K \s a. constant)) because 
in typical graphs there are not many two and four point subgraphs. 

2.6 One and Two Particle Irreducible Expansions 

Salmhofer's criterion is stated for the self-energy, i.e. the sum of all one-particle 
irreducible graphs for the two point function. Its study requires the correct 
renormalization of these contributions. Since angular sectors in a graph may 
vary from one propagator to the next in a graph, and since different sectors 
have different decays in different directions, we are in a delicate situation. In 
order to prove that renormalization indeed does the good that it is supposed to 
do, one cannot simply rely on the connectedness of these self- energy graphs, 
but one must use their particle irreducibility explicitly. So the proof requires 
a constructive particle irreducible analysis of the self-energy. The following 
the-orem summarizes the results of |541 155] : 

Theorem 2.1 The radius oj convergence of the jellium two-dimensional model 
perturbative series for any thermodynamic function is at least c/| logTj, where 
T is the temperature and c some numerical constant. As T and A jointly tend to 
in this domain, the self-energy and its first two momentum derivatives remain 
uniformly bounded so that the model is a Fermi liquid in the sense of Salmhofer. 

In the case of the jellium model J2, this analysis can be performed at the 
level of one-particle irreducible graphs [S^- The half-filled Hubbard model, 
however, is more difficult. Although there is no real divergence of the self-energy 
(the associated counterterm is zero thanks to the particle hole symmetry of the 
model at half-filling) one really needs a two-particle and one-vertex irreducible 
constructive analysis to establish the necessary constructive bounds on the self- 
energy and its derivatives [57]. For parity reasons, the self-energy graphs of the 
model are in fact not only one-particle irreducible but also two particle and one 
vertex irreducible, so that this analysis is possible. 

This analysis leads to the explicit construction of three line-disjoint paths for 
every self-energy contribution, in a way compatible with constructive bounds. 
On top of that analysis, another one which is scale-dependent is performed: after 
reduction of some maximal subsets provided by the scale analysis, two vertex- 
disjoint paths are selected in every self-energy contribution. This construction 
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allows to improve the power counting for two point subgraphs, exploiting the 
particle-hole symmetry of the theory at half-filling, and leads to the desired 
analyticity result. 

Finally an upper bound for the self energy second derivative is combined 
with a lower bound for the explicit leading self energy Feynman graph |58j . 
This completes the proof that the Hubbard model violates Salmhofer's criterion, 
hence is not a Fermi liquid, in contrast with the jellium two dimensional model. 
More precisely the following theorem summarizes the results of [56l \57\ [58] 

Theorem 2.2 The radius of convergence of the Hubbard model perturbative se- 
ries at half-filling is at least c/\ logTj , where T is the temperature and c some 
numerical constant. As T and A jointly tend to in this domain, the self-energy 
of the model does not display the properties of a Fermi liquid in the sense of 
Salmhofer, since the second derivative is not uniformly bounded. 

We would like now to enter into more technical detail, without drowning 
the reader. Hence we shall limit ourselves here to the non-perturbative analysis 
of connected functions in a single RG scale, which is the core mathematical 
problem. We compare the various models J2, J3 and H2 to the toy model of 
section |1.16( and explain why sector analysis plus momentum conservation is 
suited to analyze the two-dimensional models but fails in three dimensions. 



We claim that the J2 model in a slice is roughly similar to the Toy Model, with 
dimension d — 3, provided the momentum slice is divided into angular subpieces 
called sectors, which play the role of colors. 

The naive estimate on the slice propagator is (using integration by parts) 



(using Gevrey cutoffs fj to get fractional exponential decay). The prefactor 
corresponds to the volume of integration of the slice, M~^^ , divided by 
the slice estimate of the denominator, M"-' . This is much worse than the factor 
ji^-3j/2 ^-^g^^ would be needed. 

But the situation improves if we cut the Fermi slice into smaller pieces (called 
sectors). Suppose we divide the j-th slice into sectors, each of size roughly 
in all three directions. 

A sector propagator C-''" has now prefactor M^"^^ corresponding to the vol- 
ume of integration of the sector M^^^ divided by the slice estimate of the 
denominator, M^^ . Using integration by parts and Gevrey cutoffs fja for frac- 
tional power decay we get without too much effort the bound 



2.7 2D Jellium Model: Why Sectors Work 



(2.4) 



\C,Mx,y)\ < 6abM-^'e-^^'''^--yn 



1/2 



(2.5) 



But since N 



(2.6) 
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so that this bound is identical to that of the toy model. 

It remains just to explain why the interaction of the model is approximatley 
of the vector type. This is because of the momentum conservation rule at every 
vertex. 

In two dimensions a rhombus (i.e; a closed quadrilateral whose four sides 
have equal lengths) is a parallelogram. Hence an approximate rhombus should 
be an approximate parallelogram. 

Momentum conservation S{pi + P2 + Pa + Pi) at each vertex follows from 
translation invariance of J2. Hence Pi,P2,P3,P4: form a quadrilateral. For j 
large we have \pk \ — •\/2M/i, hence the quadrilateral is an approximate rhombus. 
Hence the four sectors to which pi, p2, Pa and p4 should be roughly equal two 
by two (parallelogram condition). 

It means that the interaction is roughly of the color (or Gross-Neveu) type 
with respect to these angular sectors: 

{J2Ma){J2Mb) (2.7) 

a b 

In fact this "rhombus rule" is not fully correct for almost degenerate rhom- 
buses. The correct statement is 

Lemma 2.3 Fix m e Z'^ . The number of A-tuples {Si, ■ ■ ■ S'4} of sectors for 
which there exist /c^ € M^, i = 1, ■ • • , 4 satisfying 

k[eS,, |fc, - fc^l < const M-J' , i = •••,4 (2.8) 

and 

\ki+ ■■■ + fc4| < const (1 + ImDM"-'' (2.9) 

is bounded by 

const(l + \m\)^M^^ {1 + j} . (2.10) 

The 1 + j factor is special to dimension 2 and is the source of painful tech- 
nical complications which were developed by Feldman, Magnen, Trubowitz and 
myself. 

The solution uses in fact anisotropic angular sectors, which are longer 

in the tangential direction (of length M^^^^). The corresponding propagators 
still have dual spatial decay because the sectors are still aprroximately flat. 

Ultimately the conclusion is unchanged: the radius of convergence of J2 in 
a slice is independent of the slice index j. 

2.8 Why Sectors Fail in = 3 

In 3D sectors and Gram's bound fail by a full power per vertex! 

There is indeed no rhombus rule in d = 3. A closed quadrilateral with equal 
sides is not a parallelogram because it can be non-planar; hence it is obtained 
by rotating half of a planar parallelogram around the diagonal by an arbitrary 
twisting angle. Therefore the jellium model interaction is not of the vector type. 
More precisely the analog of Lemma |2.3| is, with similar notations 
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Lemma 2.4 The number of A-tuples {5*1, • • • S'4} of sectors for which there 
exist g M^, i — 1, • • • , 4 satisfying 

k'^eS^, |fc, - fc^l < const M~^' , i = l, •••,4 (2.11) 

and 

\ki+ ■■■ +fc4| < const(l + |m|)M--'' (2.12) 

is bounded by 

const(l + ImD^M^J'. (2.13) 



Remark the absence of the log factor that was present in ci = 2. But for d = 3 we 
find M^^ 4-tuples which correspond to the choice of two sectors (M^-' x ) 
and of one angular twist AP . The power counting corresponds to M~^^ per 
sector propagator. Two propagators pay for one vertex integration (M^-') and 
one sector choice (M^^ ) but there is nothing to pay for the angular twist. Going 
to anisotropic sectors is possible but there remains still in this case a /2 twist 
factor. After many years of effort, we concluded that the sector method and 
Gram bound apparently cannot be improved to do better. Hence although J3 
is expected to be a Fermi liquid in the sense of Salmhofer, new methods have 
to be developed to treat it constructively [551 EZl EH] • 



2.9 The Hadamard method in a;-space 

The idea of using Hadamard's inequality in x-space to overcome the constructive 
power counting problem of Fermions in three dimensions took Jacques Magnen 
and myself four years of continuous hard work, with dozens and dozens of various 
failed trials, from 1991 to 1995 [55]. Another four or five years took place to 
fine-tune this idea for the multidimensional case with M. Disertori [67] . Another 
ten years have passed to find a momentum-conserving version, which should at 
last allow for the proof of Salmhofer's criterion in the three dimensional jellium 
model [55]. 

Let us describe the main idea in an informal way. 



We remember that the naive estimate (2.4) is far from sufficient for correct 
power counting. But we also know from perturbative power counting in momen- 
tum space that the theory should be just renormalizable. Therefore the j-slice 
propagator should behave more as the one of an infrared theory, hence 
should be bounded by 

KM-^^e-^^'''\-~y\^"'\ (2.14) 

In fact this is not totally correct, because it can be shown that the j'-th slice 
propagator at almost coinciding points, hence at |x — y| ~ is not bounded by 



M "^'K Still (2.14) is correct at typical distances for the J3 propagator in the 



j-th slice. We know from eg (2.4) that these typical space-time distances should 
be |x - 2/1 ~ M^. 
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Indeed integrating over angles on the Fermi sphere leads to an additional 
l/\x — y\ decay, because 

/ sm9ded(j)e'''°''^^''-y^ ^sm\x~y\/\x-y\ (2.15) 
Jo 

Hence for typical distances ja:; — y| ~ the propagator indeed obeys the 
improved estimate 

|C,(x,y)|,_,|^M.| < i^M-2■'"e-[*^"^l--^ll''^ (2.16) 

The problem is that this bound is wrong at small distances. However at 
small distances there is a bonus, namely the corresponding integration volumes 



are smaller. Another related problem is that if we use the Gram bound (1.88 1 



to bound a determinant with all the matrix elements Cj{xp — yq) corresponding 



to large distances \xp — yq\, we still loose the improvement (2.161, because the 
norms in (1.88) will correspond again to propagators at coinciding points! 



The solution is to use the Hadamard bounds ( 1.89 )-( 1.90 1 because they con- 
serve the decay of the typical propagators. 



Remember however, as seen conveniently in (1.91) that Hadamard bounds 
consume the 1/n! symmetry factor for n vertices. Hence we can no longer use 
the explicit tree formula and the method of |46l [47] . 

But we can still use the standard old-fashioned cluster expansion between 
cubes. Summarizing: 

• Just renormalizable power counting is recovered for the main part of the 
theory if we use Hadamard's bound rather than Gram's bound. 

• The factor n\ is lost in the Hadamard bound at order n; this forces us to 
rely on the on-canonical tool of cluster expansion between cubes 

Still, this solves the constructive problem only for the main part of the 
propagator, the one at typical distances |a; — yj ~ . However at smaller 
distances there is a bonus, namely the volume factors for spatial integration are 
also smaller. 

It turns out that the problem of smaller than typical distances can be solved 
with an auxiliary superrenormalizable decomposition 



C^=^C^'= (2.17) 



fe=0 

of the propagator. Roughly speaking C^'^ corresponds to \x — y\ ~ , It means 
that even the single slice theory in d > 3 is a non-trivial theory that contains a 
rather non-trivial renormalization, like the one of (/)|. This renormalization can 
be analyzed by means of the auxiliary scales. The solution is in fact more com- 
plicated that what we sketch here, and has to take into account the anisotropy 
between the space and the imaginary time variables [66| I67j . 
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The use of non-canonical lattices of cubes in this method is the signal that we 
have probably still not found the optimal constructive treatment of J3 . Recently 
we found a better decomposition that should allow to check Salmhofer's criterion 
for J3 [55] • However in its present stage it will still use a non-canonical cluster 
expansion between cubes. It would be interesting to find a solution such as the 
loop vertex expansion that solves the constructive problem in a truly canonical 
way. 



2.10 2D Hubbard Model 

The Hubbard model lives on the square lattice 'E? , so that the three dimensional 
vector X — {xo,x) is such that x — (711,712) € From now on we write vi 
and V2 for the two components of a vector v along the two axis of the lattice. 

At half-filling and finite temperature T, the Fourier transform of the prop- 
agator of the Hubbard model is: 

Cab{k)^Sab-, '^—TTT, e(k) = COS /ci -f COS fc2 , (2.18) 

iko - e[k) 

where a, 6 e {'f , J,} are the spin indices. The vector k lives on the two-dimensional 
torus IR^/(27rS)^. Hence the real space propagator is 

Cabix) = j / e*'" Cafc(fc) . (2.19) 

Recall that |fco| > i^/P 7^ hence the denominator in C{k) again can never 
be at non zero temperature. This is why the temperature provides a natural 
infrared cut-off. When T — ^ (which means /3 — >■ cx)) fcg becomes a continuous 
variable, the discrete sum becomes an integral, and the corresponding prop- 
agator Cq(x) becomes singular on the Fermi surface defined by A:o = and 
e(k) = 0. This Fermi surface is a square of side size v^tt (in the first Brillouin 
zone) joining the corners (±7r, 0), (0, ±7r). We call this square the Fermi square, 
its corners and faces are called the Fermi faces and corners. Considering the 
periodic boundary conditions, there are really four Fermi faces, but only two 
Fermi corners. 

In the following to simplify notations we will write: 

I (fk = j d?k , I d?x = ]- r dxo . (2.20) 

In determining the spatial decay we recall that by anti-periodicity 

Cix) = /(xo,x) :^ Y (-1)" {xo + ^,x) . (2.21) 

where Co is the propagator at T = 0. Indeed the function / is anti-periodic and 
its Fourier transform is the right one. 
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The interaction of the Hubbard model is again (2.3 1: 

Sv = x f d^x (Vv^V)'(^) , 

Jv 



(2.22) 



where V :— [—f3, f3] x V' and V is an auxiliary finite volume cutoff in two 
dimensional space that will be sent later to infinity. 



2.11 Scale Analysis 

The theory has a natural lattice spatial cutoff. To implement the renormal- 
ization group analysis, we introduce as usually a compact support function 
u{r) G Cg°(R) (it is convenient to choose it to be Gevrey of order a < 1 so as 
to ensure fractional exponential decrease in the dual space) which satisfies: 

u{r) = for |r| > 2 ; u{r) = 1 for |r| < 1 . (2.23) 

With this function, given a constant M > 2, we can construct a partition of 
unity 

oo 

1 = ^?i.(r) Vr^O ; 

uo{r) = 1 - u{r) ; u,(r) = u{M^''''~^^r) - u{M^W) for z > 1 . (2.24) 
The propagator is then divided into slices according to this partition 

oo 

C(fc)-^a(fc) (2.25) 

i=0 

where 

a{k) = C{k)u,[kl + e^{k)] . (2.26) 

(indeed ki + e^{k) > > 0). 

In a slice of index i the cutoffs ensure that the size of fcp + (k) is roughly 
M~^'. More precisely in the slice i we must have 

M"2i < fcg + e^{k) < 2M^M''^' . (2.27) 

The corresponding domain is a three dimensional volume whose section 
through the fco = plane is the shaded region pictured in Figure [6j 

Remark that at finite temperature, the propagator Ci vanishes for i > 
imaxiT) where M*— ~ 1/T (more precisely i^a^T) = £;(log ^/logAf), 
where E is the integer part), so there is only a finite number of steps in the 
renormalization group analysis. 

Let us state first a simple result, for a theory whose propagator is only Ci, 
hence corresponding to a generic step of the renormalization group: 
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Figure 6: A single slice of the renormalization group 



Theorem 1 The Schwinger functions of the theory with propagator Ci and 
interaction (2.22) are analytic in \ in a disk of radius Ri which is at least c/i 



for a suitable constant c: 



Ri > c/i 



(2.28) 



The rest of this section is devoted to the definitions and properties of H2 
sectors, their scaled decay and momentum conservation rules. As discussed 
already this result is a first step towards the rigorous proof [571 EH] that H2 is 
not a Fermi liquid in the sense of Salmhofer. 



2.12 Sectors 

This section is extracted from [56] . 

The "angular" analysis is completely different from the jellium case. We 
remark first that in a slice, kg + e^{k) is of order M~^*, but this does not fix 
the size of e^(k) itself, which can be of order M~^^ for some j > i. In order 
for sectors defined in momentum space to correspond to propagators with dual 
decay in direct space, it is essential that their length in the tangential direction is 
not too big, otherwise the curvature is too strong for stationary phase methods 
to apply. This was discussed first in [62]. This leads us to study the curve 
(cos/ci + cosA:2)^ = M~'^^ for arbitrary j > i. We can by symmetry restrict 
ourselves to the region < fci < tt/2, k2 > 0. It is then easy to compute the 
curvature radius of that curve, which is 

R= i-^'k^+--^'k2f'\ . (2.29) 
I cos fci sin ^2 + cos fc2 sin fci | 

We can also compute the distance d{ki) to the critical curve cos fci +cos fc2 = 0, 
and the width w{ki) of the band M^^ < \ cos fci +cos fc2| < V^M.M^^. We can 
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then easily check that 



dih) - wik,) ^ j^^^^ , (2.30) 
Rih) ^ -^^^ — , (2.31) 

where f ^ g means that on the range Q < ki < tt /2 we have inequahties 
cf < g < df for some constants c and d. 
Defining the anisotropic length 



/(fci) = ^w{ki)R{ki) ~ + fci , (2.32) 

the condition in [52] is that the sector length should not be bigger than that 
anisotropic length. This leads to the idea that fci or an equivalent quantity 
should be sliced according to a geometric progression from 1 to M~^l'^ to form 
the angular sectors in this model. 

For symmetry reasons it is convenient to introduce a new orthogonal but 
not normal basis in momentum space (e+,e_), defined by e+ = (l/2)(7r,7r) and 
e_ = (l/2)(— vr, tt). Indeed if we call (fc+,fc_) the coordinates of a momentum 
k in this basis, the Fermi surface is given by the simple equations kj^ ~ ±1 or 
k- = ±1. This immediately follows from the identity 

cosfci +cos/c2 — 2 cos(7r/c+/2) cos(7rfc_/2) . (2.33) 

(Note however that the periodic b.c. are more complicated in that new ba- 
sis). Instead of slicing e(k) and fci, it is then more symmetric to slice directly 
cos(7rfc+/2) and cos(7rfc_/2). 

Guided by these considerations we introduce the partition of unity 

vo{r) = 1 - u{M^r) 
Vs{r) ; { vs ^ Us+i forl<s<i-l (2.34) 

Vi{r) = m(M2V) 



and define 



where 



Q(fc)- J2 C.^^k) (2.35) 
Claik) = Q{k)vs^ [cos2(7rfc+/2)] [cos^ 7rfc_/2)] . (2.36) 



We remark that using (2.27) in order for Ci^a- not to be 0, we need to have 
s+ + s_ > i — 2. We define the "depth" l{a) of a sector to be Z = s+ +s_ — i + 2. 

To get a better intuitive picture of the sectors, we remark that they can be 
classified into different categories: 

• the sectors (0,i) and (i,0) are called the middle-face sectors 

• the sectors (s,i) and (i,s) with < s < i are called the face sectors 
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• the sector (i,i) is called the comer sector 

• the sectors (s,s) with (i — 2)/2 < s < i are called the diagonal sectors 

• the others are the general sectors 

Finally the general or diagonal sectors of depth for which s+ + s_ = i — 2 
are called border sectors. 

If we consider the projection onto the {k-^-,k-) plane, taking into account 
the periodic b.c. of the Brillouin zone, the general and diagonal sectors have 
8 connected components, the face sectors have 4 connected components, the 
middle face sectors and the corner sector have 2 connected components. In the 
three dimensional space-time, if we neglect the discretization of the Matsubara 
frequencies, these numbers would double except for the border sectors. 

2.13 Scaled decay 

Lemma 2.5 Using Gevrey cutoffs of degree a < 1, the propagator Ci_a obeys 
the scaled decay 

|C,,„| < c.Af-*^'e-[''''"(^'y)l° (2.37) 

where 

d^M^, y) = {M~'\xo - yol + A/"^+ \x+ -^y+\+ M"'- |x_ - y-\} . (2.38) 



Proof This is essentially Fourier analysis and integration by parts. If a; = 
(ni, 712) S we define (a;_|_, a;_) — (7r/2)(7ii+n2, —"-i)- The vector (x_|_, a:_) 
then belongs to (7r/2)S^ but with the additional condition that X-^. and X- have 
the same parity. 



Defining, for X e [{n/2)W, f 

A,.(^)= (i/2)^E/7'^^+/ 



+2 /.+2 

dk g*(*^o2;o+fe+a;+ + fe_x_ ) 



-2 

U2 



Ui[k^ +4cos2(7rfc+/2)cos2(7rfc_/2)] 
iko — 2cos(7rA:+/2) cos(7rfc_/2) 

^ [cos2(7rfc+/2)] [cos2(7rfc_/2)] (2.39) 



V, 



we note that Ci^a{X) = Di^„{X) for X satisfying the parity condition. 

(Remember the Jacobian ^ from dkidk2 to dk^dk^, and the initial domain 
of integration that is doubled.) 

The volume of integration trivially gives a factor M~* for the fco sum and 



factors M~^+ and M~^~ for the fc+ and k- integration (see (2.44) below). 
The integrand is trivially bounde d by on the integration domain, and this 
explains the prcfactor cM^^^^ in (2.371. 
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We then apply standard integration by parts techniques to formulate the 
decay. From e.g. Lemma 10 in [51] we know that to obtain the scaled decay of 
Lemma |2.5| we have only to check the usual derivative bounds in Fourier space: 



:A,all < A.B"iVr"»M''+"+M''-"-(n!)i/" (2.40) 



"9fco« dkl+ dkg 

where n — hq + n+ + n_, and the derivative really means the natural 
finite difference operator (l/27rT)(/(fco + 27rT) — /(fco)) acting on the discrete 
Matsubara frequencies. The norm is the ordinary sup norm. 



But from ( |2.39[ ), 

^ _ 1 u^/cg +4cos2(7rfc+/2)cos2(7rfc_/2)] 

''"^ ' ~ 16/3 iko~2 cos(7rfc+/2) cos(7rfc_/2) 

?;s^[cos^(7rA:+/2)] t)s_ [cos^(7rfc_/2)] (2.41) 

and the derivatives are bounded easily using the standard rules for derivation, 
product and composition of Gevrey functions, or by hand, using the support 
properties of the v^^ and Vs_ fonctions. For instance a derivative can act on 
the Vs^ [cos^(7rfc+/2)] factor, in which case it is easily directly bounded by cM^+ 
for some constant c. When it acts on u^^o + 4cos^(7rfc+/2) cos^(7rA;_/2)] it is 
easily bounded by c.M^*~*'+~^''- hence by c.M^+ , using the relation s+ + s_ > 
i — 2. When it acts on the denominator [iko — 2 cos(7rfc_|_/2) cos(7rA;_/2)]^^. it is 
bounded by c.M*~''~ , hence again by c.M''+ , using the relation s+ + s_ > i — 2. 
Finally when it acts on a cos(7rA;+/2) created by previo us derivations, it costs 
directly c.M'+ . The factorial factor (?i!)i/" in (2.40) comes naturally from 



deriving the cutoffs, which are Gevrey functions of order a; deriving other factors 
give smaller factorials (with power 1 instead of 1 /a) . 

Finally a last remark: to obtain the Lemma for the last slice, i = imax{T), 
one has to take into account the fact that xq lies in a compact circle, so that 
there is really no long-distance decay to prove. 

2.14 Support Properties 

If Ci^a{k) 7^ 0, the momentum k must obey the following bounds: 

|/co| < (2.42) 

'M-1 < I cos(7rfc±/2)| < 1 for s± = , 

M-^±-i < I cos(7rfc±/2)| < \/2M-^± for 1 < s± < i - 1 , (2.43) 
I cos(7rfc±/2)| < ^/2M-'' for s± = i . 

In the support of our slice in the first Brillouin zone we have |fc+| < 2 and 
< 2 (this is not essential but the inequalities are strict because i > 1)- It 
is convenient to associate to any such component k± a kind of "fractional part" 



47 



called q± defined hy q± ^ k± — 1 ii k± > and q± = k± 
^ \q± \ < 1- Then the bounds translate into 



1 if fc± < 0, so that 



'2/7rM < |g±| < 1 
2M-'±/nM <\q±\< V^M' 
\q±\ < 



for s± = , 

for 1 < s± < j — 1 

for s± = i . 



(2.44) 



2.15 Momentum conservation rules at a vertex 

Let us consider that the four momenta fci, ^2, k^, k^, arriving at a given vertex 
V belong to the support of the four sectors cri, (T2, fs, a^, in slices ii, 12, is, ?4. 



In Fourier space the vertex (2.221 implies constraints on the momenta. Each 



spatial component of the sum of the four momenta must be an integer multiple 
of 27r in the initial basis, and the sum of the four Matsubara frequencies must 
also be zero. 

In our tilted basis (e-|-,e_), this translates into the conditions: 



ki.o + k2fi + /cs^o + fc4,o = , 

+ fc2,+ + fc3^+ + /S4,+ = 2n4 

ki^- + ^2,- + fca.- + ki^- — 2n_ 



(2.45) 
(2.46) 
(2.47) 



where and n_ must have identical parity. 

We want to rewrite the two last equations in terms of the fractional parts 
91, 92, 93 and 54- 



Since an even sum of integers which are ±1 is even, we find that (2.46) and 



(2.471 imply 



9i, 
91, 



92. 
92, 



93,H 

93.- 



94,+ = 2rn+ 
(74__ = 2m_ 



(2.48) 
(2.49) 



with m+ and m_ integers. Let us prove now that except in very special cases, 
these integers must be 0. Since |9j.±| < 1, \m±\ < 2. But \qj,± \ = 1 is possible 
only for Sj^± — 0. Therefore \m±\ = 2 implies Sj_± = Vj. Now suppose e.g. 
|to+| = 1. Then Sj^+ is for at least two values of j. Indeed for sj _+ 7^ we 
have |9j,±| < V2M~^, and assuming 3-\/2M^^ < 1, equation (2.48) could not 
hold. 

We have therefore proved 



unless Sj _|_ is for at least two values of j , and rn_ 



Lemma 2.6 

unless Sj^- is for at least two values of j. 







Let us analyze in more detail equations (2.48) and (2.49) for |m+| = |m_| = 
0. Consider e.g. (2.48). By a relabeling we can assume without loss of generality 

< ii, in which case 



that si_+ < S2,+ < S3,+ < S4,+ Then either si_+ — ii or si, 
combining equations (2.48) and ( |2.44 1 we must have: 



3y2M"''^'+ > 2Af"''i'+/7rM 



(2.50) 
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which means 



S2,+ < + 1 + 



log(37r/V2) 
logM 



(2.51) 




(2.52) 



The conclusion is: 

Lemma 2.7 If m± ~ 0, either the smallest index Si^± coincides with its scale 
ii, or the two smallest indices among Sj ± differ by at most one unit. 

Now we can summarize the content of both Lemmas in a sUghtly weaker but 
simpler lemma: 

Lemma 2.8 

A) (single slice case) 

The two smallest indices among Sj,+ for j — 1,2,3,4 differ by at most one 
unit, and the two smallest indices among Sj „ for j = 1,2,3,4 differ by at most 
one unit. 

B) Multislice case 

The two smallest indices among + for j — 1,2,3,4 differ by at most one 
unit or the smallest one, say si^+ must coincide with its scale ii, which must 
then be strictly smaller than the three other scales i2, ia and i^. Exactly the 
same statement holds independently for the minus direction. 

This lemma allows to check that again the single scale analysis works and 
leads to a radius of convergence independent of the slice [S^ . 

2.16 Multiscale Analysis 

The half-filling point is very convenient since the particle-hole exact symmetry 
at this point ensures that there is no flow for the Fermi surface itself. 

Using the sector decomposition and the momentum conservation we find 
that power counting of the 2D Hubbard model is essentially similar to one 
dimensional case with logarithmic corrections |57j . Expanding the two point 
function to second order we can even find lower bounds which prove that this 
model is not a Fermi liquid in the sense of Salmhofer [SS] . 
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